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Abstract 


Modeling of Thermal Lensing in Side and End-Pumped 
Finite Solid-State Laser Rods 

Student: Vincent G. Brackett Department: Mathematics 

Date of Birth: December 17, 1942 Advisor: Dr. U.O.Farrukh 

An analytical expression for approximating the time- 
dependent thermal focal length in finite solid-state laser 
rods has been derived. The analysis is based on the tem- 
perature variation of the material refractive index caused 
by optical pumping of these rods. Several quantities were 
found to be relevant to this analysis. These quantities 
were the specific thermal profiles of the rods, type of op- 
tical pumping employed, type of cooling scheme employed 
(side and end-cooling parameters), and the specific materi- 
al characteristics of the rods. The Thermal Lensing Model 
was formulated using the geometric ray tracing approach. 

The focal lengths are then approximated, by calculating the 
phase shift in the index of refraction, as the different 
rays of an incident plane wave are tracked through a lens- 
like crystal medium. The approach also applies in the case 
of Gaussian or parabolic pump beams. It is shown that the 
prediction of thermal focal length is in good quantitative 
agreement with experimentally obtained data. 
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CHAPTER 1 


Introduction 

The next generation of solid-state lasers (to be used 
for remote sensing of the atmosphere from space platforms), 
are under development. Various new lasing materials and 
laser conf iguati ons are being investigated. Currently, the 
solid-state, end-pumped diode lasers are the most promising 
candidates due to their small size, reliability, perform- 
ance, low power and cooling requirements, etc. Analyzing 
thermal lensing effects for these end-pumped conf iguations , 
using various cooling schemes are of interest to the de- 
signers of laser systems. Theoretical formulation, numeri- 
cal calculations, and experimental data on thermal lensing 
does not exist for end-pumped conf iguations . 

The management of thermal energy deposited in a solid- 
state material is critical in the design of laser systems. 1 
Any solid laser material operating in either pulsed or con- 
tinuous wave (CW) mode of operation dissipates an apprecia- 
ble amount of heat. The heat arises from non-radiative 
energy transitions as the pump beam traverses the material . 
Stated very briefly, when a material is excited to provide 
more atoms (or molecules) in a higher energy state than in 
some lower state, the material will be capable of ampli- 
fying radiation at the frequency corresponding to the ener- 
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gy level difference . 2 The acronym LASER derives its name 
from this process: "Light Amplification by Stimulated Emis- 
sion of Radiation." 

It’s during this process that some fraction of the to- 
tal energy of the pump beam is absorbed by the lasing mate- 
rial. This absorbed energy is then dissipated as heat. In 
the cylindrical geometries normally used, the heat is ex- 
tracted from the surface of the cylinder, thereby gener- 
ating a radial thermal gradient. Both radial and axial 
gradients are present in end-pumped lasers. Employment of 
active cooling to dissipate heat from laser crystals will 
also affect the thermal gradients. Therefore, the thermal 
gradients present in laser rods are governed by many fac- 
tors. These include: geometries of the rod in use, type of 
laser material, thermal diffusivity of the material, and 
type of pumping and active cooling schemes employed. A 
Temperature Distribution Model has been developed to model 
the time dependent temperature distributions of finite 
solid-state laser rods that include the above factors. ^'3 
The change in temperature, within a laser rod, causes 
thermal defocusing and distortion of the laser beam. This 
is due to the temperature variation of the material refrac- 
tive index. These thermal gradients cause stresses within 
the laser material. The generated stresses induce bire- 
fringence and non-uniform expansion (bulge) in the laser 
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rod. Collectively, the above factors of temperature varia- 
tion of the refractive index, induced birefringence, and 
bulge are called thermal lensing. Thermal lensing is the 
defocusing and distortion of a laser beam transmitted 
through a material medium, due to the non-uniform heating 
of the material. ^ 

In this thesis, using the Temperature Distribution 
Model, thermal lensing effects will be analyzed to develop 
a Thermal Lensing Model. The Thermal Lensing Model will be 
employed to predict changes in the focal length and resul- 
tant distortions due to the time dependent temperature var- 
iations of the refractive index for two common laser ma- 
terials that are under evaluation. These materials are 
Neodymium doped Yttrium Aluminum Garnet (YAG) (Y3AL5O12) 
and Titanium doped Aluminum Oxide (Sapphire) (AL2O3). From 
this point onward, these two materials will be called 
Nd:YAG and Ti:AL203 respectively. They will be evaluated 
for various pumping (side and end) and cooling schemes. 

The induced birefringence and bulge caused by stress in the 
laser material will be discussed. However, these factors 
will not be implemented in the first version of the Thermal 
Lensing Model as they are minor in nature. 

In the next chapter, the Temperature Distribution Mod- 
el will be discussed. It's from this model's output that 
the time dependent temperature profiles are derived. These 
profiles serve as input to the Thermal Lensing Model. 
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CHAPTER 2 


Temperature Distribution Model 

The problem of managing the thermal energy deposited 
within a solid-state laser material under active pumping 
has long been of critical importance in the design of laser 
systems. 1 Koechner's paper on the transient temperature 
profile in solid-state lasers was the original work in this 
area. 5 i n it, the time dependent temperature distribution 
of an infinitely long laser rod with a uniform heat source 
was determined. This configuation was adequate to treat 
the flashlamp pumping of relatively long laser rods, and 
rods with little or no end cooling employed. However, this 
treatment has become inapplicable to other pump configua- 
tions. Recently, several new solid-state laser materials 
have been introduced, each with their own thermal deposi- 
tion mechanisms . ® Additionally, several other novel laser 
techniques such as end-pumping with laser diodes have also 
been developed. For these systems, the assumptions of a 
uniform distribution of pump energy and an infinite rod 
length are no longer valid. The one dimensional approach 
of Koechner 5 does not apply when the length of the laser 
rod is comparable to its diameter. Furthermore, the longi- 
tudinal symmetry of an infinite rod is broken when it is 
end-pumped . 
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It was for these reasons that the Temperature Distri- 
bution Model was developed. The time dependent temperature 
distribution of a cylindrical rod of finite length has been 
derived by Farrukh et al . in two separate papers. -*-'3 Addi- 
tionally, they studied the evolution of thermal profiles 
for various pumping and cooling schemes. The general ex- 
pression of the temperature distribution is given as:' 

oo 

T(r.z) = ^ * f W (2-1) 

n,m=l 

where 0 n is the radial eigenfunction , 'y , m is the axial ei- 
genfunction, C njn is the pump pulse expansion coefficient, 
and f(t) is the time dependence function. The time depen- 
dence function f(t) for the on and off periods of the 
pumping pulse is given, respectively as: 7 
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where K is the inverse of diffusivity, Lamda is the thermal 
conductivity, and x nm is the time constant for specific 
eigenfunctions . 

The expression for temperature distribution is general 
and applies to various situations such as pulsed or contin- 
uous operations of the laser. A variety of different pum- 
ping and cooling schemes including end-pumping by a Gaus- 
sian beam were also discussed. To gain an appreciation for 
the capability of this model, we will look at two of the 
more interesting cases that this model can be used to eval- 
uate. These are end-pumped single and multiple pulse oper- 
ations for two common solid-state laser materials (Nd:YAG 
and Ti : AL2O3 ) . 

A. Single Pulse Operations 

If a laser rod is pumped by a single pump pulse, a 
transient thermal profile will be established. This is the 
result of a fast pump-introduced heating process followed 
by a slow recovery of the rod to its ambient temperature. 
Although this thesis will primarily address long pump pul- 
ses, it is appropriate to distinguish the difference be- 
tween short and long pump pulses. If the duration of the 
pumping pulse is very short in comparison to the laser rod 
relaxation time, then the heat energy deposited in the rod 
is considered to be deposited instantaneously. (The laser 
rod relaxation time is defined as the time required for the 
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rod to relax to its ambient temperature.) Rod relaxation 
times are approximately 3.5 seconds for Nd:YAG and 300ms 
for Ti : AL 2 O 3 for typical rod dimensions. When a pulse is 
characterized as short or long, this characterization de- 
pends on the laser rod material. If the pump pulse dura- 
tion is less than one percent of the relaxation time of the 
rod, then the pump pulse is considered to be short. In 
this case, as previously stated, the heat energy is depos- 
ited instantaneously. An example of this will be shown in 
a later graph. For long pulse operations, the heat energy 
being deposited builds up and relaxes simultaneously during 
the pumping pulse. 

Figure 2-1 is an illustration of the transient nature 
o.f temperature distribution in an optically pumped cylin- 
drical laser rod for a long pumping pulse. During the 



Figure 2-1. Thermal Relaxation of a Single Input 
Pulse in a Cylindrical Laser Rod^ 


7 



pulse, heat energy is absorbed by the laser crystal. At 
the end of the pump pulse, the temperature of the rod has 
increased above the ambient in a spatially non-uniform man- 
ner. Since only the surface of the rod is in contact with 
the coolant, a radial and axial heat flow develops during 
and at the end of the pump pulse. The initial non-uniform 
spatial temperature distribution changes to a nearly para- 
bolic temperature profile which then decays at a rate de- 
termined by the rod specific relaxation time constant. 

Since lasing action occurs during or at the end of the pump 
pulse, the optical distortion in single pulsed lasers 
arises from the thermal gradients established during the 
pumped cycle. 

Figure 2-2 and 2-3 are examples of the Temperature 
Distribution Model output for the time dependent heat prop- 
agation in end-pumped Nd : YAG and Ti:AL 2 0 3 laser rods for a 
pump beam with spatially circular and Gaussian profiles re- 
spectively, with moderate side cooling. In both of these 
figures, the energy absorbed as heat, is lJ/cm 3 at the 
front end of the rod for the circular pump and at the cen- 
ter of the front end for the Gaussian pump. Figure 2-2 was 
chosen as an example because it shows the impact of side 
cooling in generating radial gradients in a Nd : YAG rod af- 
ter only 20ms into its thermal relaxation cycle. Also in 
Figure 2-2, the temperature at the center of rod has not 
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appreciably relaxed toward the ambient temperature at the 
20ms point in time. Figure 2-3 was chosen as an example to 
illustrate the difference between a circular and Gaussian 
pumping beam. This figure also shows that after only 20ms 
into the thermal relaxation cycle, the Ti:AL 203 laser rod 
has significantly relaxed toward its ambient temperature. 
However, the thermal energy content in the Gaussian pump 
case is only about 1/8 of that of the circular pump case. 

Figure 2-4 is an example of the model output of the 
temperature rise above ambient for five different pulse 
durations for a Nd:YAG laser rod with one pulse of 50mJ of 
absorbed energy. The Nd:YAG laser was end-pumped with a 
Gaussian beam and subjected to moderate side cooling. The 
rod ends were not cooled. The 1ms pulse in this figure is 
an example of a short pulse where the heat energy is depos- 
ited instantaneously. The remaining pulses in this figure 
are examples of relatively long pulses. In Figure 2-4, the 
rod has essentially relaxed to the ambient temperature ap- 
proximately 3.5 seconds after completion of the pumping 
pulse. 

Figure 2-5 is an example of the model output of the 
temperature rise above ambient for five different pulse 
durations for a Ti:AL 2 C >3 laser rod using 50mJ of absorbed 
energy per pulse. The Ti:AL 2 C >3 laser rod was end-pumped 
with a Gaussian beam and subjected to moderate side- 


11 



TEMP CHANGE (deg C) TEMP CHANGE (deg C) 


4.00 


Single Pulse 

Energy =* 50rrW (absorbed) 



Figure 2-4. Temperature Profiles for Various 

Durations of a Single Pump Pulse for 
Nd:YAG (At Center of Pump-End Face) 



Figure 2-5. Temperature Profiles for Various 

Durations of a Single Pump Pulse for 
TitALoOj (At Center of Pump-End Face) 
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cooling. The rod ends were not cooled. The 1ms pulse in 
this figure is an example of a short pulse where the heat 
energy is deposited instantaneously. The remaining pulses 
in this figure are examples of relatively long pulses. In 
Figure 2-5, the rod has essentially relaxed to thermal 
equilibrium approximately 300ms after completion of the 
pumping pulse. 

Comparison of Figures 2-4 and 2-5 also shows that the 
Ti:AL 203 rod has a greater temperature rise than the Nd : YAG 
rod under the same conditions of absorbed energy, pulse 
duration, and cooling scheme. These results are not unex- 
pected as the thermal conductivity of Ti:AL 2 C >3 is higher 
than Nd : YAG . 

B. Multiple Pulse Operations 

In the previous section, we discussed pump induced 
thermal profiles that occur for single-pulse operations. 

The pulse input intervals are long with respect to the 
thermal relaxation times of the laser rod. Additionally, 
the rod has returned to the ambient temperature before the 
onset of the next pumping pulse. If the laser is operated 
at a repetition rate at which there exists a residual tem- 
perature distribution from the previous pulse, a tempera- 
ture buildup in the rod will occur. The temperature at 
specific points will increase until a steady-state condi- 
tion is reached. Figure 2-6 is an illustration of the tem- 
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Figure 2-6. Thermal Relaxation in a Solid Cylinder 
with Pulse Input Intervals near the 
Thermal Relaxation Times^ 


perature profile of a laser operated at a low repetition 
rate . 

Figure 2-7 is an illustration of thermal buildup at a 
pulse interval smaller than the thermal relaxation time. 
Under these conditions, the laser will soon establish a 
steady-state operating condition where the radial and axial 
rate of flow of heat out of the rod will equal the rate of 
heat deposition. During the transition period, the rod 
radial and axial thermal profiles will continue to' change 
until the steady-state equilibrium temperature is reached. 
In the extreme case, when the interval between pulses be~ 
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Figure 2-7. Thermal Buildup when Pulse Interval — is_ 

Small Compared to Thermal Rela xation Time^ 

comes negligibly small compared to the thermal relaxation 
time, the thermal profile of the laser rod will approach 
the profiles of continuously pumped rods with the same to- 
tal power deposition. The transient temperature distribu- 
tions, as the rod approaches thermal equilibrium, has a 
large effect on the laser energy output and beam divergence 
of the laser^. These transient periods will be discussed 
in the next chapter on thermal lensing. 

Figures 2-8 and 2-9 are examples of a Gaussian end- 
pumped Nd : YAG laser being operated at a high pulse repeti- 
tion rate. Figure 2-8 depicts the detailed buildup of tern- 
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perature for two pump cycles. Figure 2 9 shows the maximum 
and minimum temperatures within each pulse cycle. These 
graphs show that the Nd : YAG laser will reach a steady-state 
temperature in approximately 3.0 seconds. The estimated 
maximum temperature of the Nd:YAG rod is approximately 
7 2 . 5°C above the ambient at the center of the rod pump-end 
face . 

This knowledge of the maximum temperature is critical 
to laser design. For example, if the temperature exceeds 
the stress fracture limit of the rod, the laser crystal 
will fracture. The stress fracture temperature limit for 
Nd : YAG is approximately 115°C. 2 

Figure 2-10 and 2-11 are examples of the output from 



Figure 2-10. Variation of Temperature for Ti:AL-?0^ 

(At f.Rnter of Pump-End Face) 


17 



o 

o> 

<D 

~0 


UJ 

CD 

X 

< 

X 

o 

CL 

LxJ 



Figure 2-11. Extreme Temperature for TirALgC^ 

(At Center of Pump-End Face) 


the Temperature Distribution Model of the temperature rise 
due to accumulation effects of multipulse operations for 
Gaussian end-pumped Ti:AL 2 0 3 rods. As can be seen, the rod 
reaches a steady-state temperature in approximately 300ms 
and has a maximum temperature of approximately 42°C at the 
center of the pump-end face. 

The above presented example profiles from the Tempera- 
ture Distribution Model are just a few of the many profiles 
that can be generated. Other profiles, such as radial and 
axial temperature profiles versus time, can also be gener- 
ated for a variety of pumping and cooling schemes. 
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C. Summary of the Temperature Distribution Model 

The general theory describing the transient thermal 
profiles in CW and pulse-pumped lasers for a large variety 
of operating parameters are presented in Koechner^/S and 
Farrukh et al.^'^/' These general theories were imple- 
mented by Farrukh into a general analytical model . The 
model then generates thermal profiles for a large variety 
of finite solid-state lasers, with user specified operating 
parameters. Examples of typical profiles were presented 
throughout this section to give the reader a feel for the 
capabilities of this model. 

These time-dependent derived thermal profiles are 
critical in determining thermal lensing effects which de- 
pend upon the spatial distribution of temperature and its 
gradients . ® 

In the next chapter, the theory of thermal lensing and 
formulation of the Thermal Lensing Model will be discussed. 
The model takes as its input, the time-dependent temper- 
ature profiles from the Temperature Distribution Model out- 
lined above. The Thermal Lensing Model is then used to 
estimate the focal lengths and distortions present in fi- 
nite solid-state Nd:YAG and Ti:AL 203 laser rods. 
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CHAPTER 3 


Thermal Lensing 

A. Introduction to Thermal Lensing. 

Thermal lensing is the distortion and defocusing of a 
laser beam transmitted through material medium, due to the 
non-uniform heating of the material. ^ Any solid-state la- 
ser material operating in either the steady-state, pulsed, 
or CW mode of operation must dissipate an appreciable 
amount of heat. This is the result of a fast pump-intro- 
duced heating process (where some fraction of the total 
pump energy is deposited in the laser rod as heat), fol- 
lowed by a slow recovery of the rod to its ambient tempera- 
ture. In the cylindrical rod geometries normally used, the 
heat is removed on the surface of the rod. This process 
generates radial and axial thermal gradients in the rod. 

The change in temperature, within a laser rod, causes 
thermal defocusing and distortion of the laser beam. This 
is due to the temperature variation of the material refrac- 
tive index. Additionally, these thermal gradients cause 
stresses within the laser material. The generated stresses 
induce birefringence and non-uniform elongation of the rod 
and bulging of the ends. Collectively, these above factors 
of temperature variation of the refractive index, birefrin- 
gence, and bulge are referred to as thermal lensing. From 
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Koechner^, the temperature dependent variation of the re- 
fractive index constitutes the major contribution to ther- 
mal lensing. The stress dependent variation of the re- 
fractive index (birefringence) modifies the focal length 
about 20 percent . ^ The effects of end-face curvature, 
caused by an elongation of the rod on thermal lensing is 
less than six percent. ^ 

The resulting distortion and defocusing of laser beams 
have been observed in a variety of experiments. Several 
theoretical treatments of thermal lensing have been given. 
These include geometrical -optics-*-® , wave-optics-*--*-, and the 
vector Kirchhoff diffraction formulation of thermal len- 
sing^'-*-! to account for thermally-induced birefringence. 

In the general theory for the formulation of the Ther- 
mal Lensing Model, various effects of thermal lensing will 
be discussed. However, in the model, only thermal lensing 
due to the temperature variation of the material refractive 
index using ray-trace methodology will be implemented. The 
ray-trace method was chosen as the simplest approach to ob- 
tain a working model . This method makes maximum use of the 
temperature profiles generated by the Temperature Distribu- 
tion Model previously discussed. This approach also ac- 
counts for 75-80 percent of the thermal lensing effects 
present in side and end-pumped laser rods. In general, 
the model will give the user a first approximation of the 
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thermal focal length and the degree of distortion expected 
for rods of specific dimensions and specified operating pa- 
rameters. The input to this approach will be the tempera- 
ture profiles generated by the Temperature Distribution 
Model . 

In the next section, a detailed one dimensional theory 
of thermal lensing from Koechner 2 ' 9 will be presented. The 
presented theory will outline the effects of temperature, 
birefringence, and bulge that contribute to thermal lensing 
for a side-pumped Nd:YAG laser rod. 

B. One Dimensional Theory of Thermal Lensing. 

The one dimensional theory of thermal lensing is 
looked at because the basic underlying theory and mathemat- 
ical formulation can be modified and applied to the problem 
of developing the Thermal Lensing Model. 

Consider the case of the CW side-pumped Nd:YAG rod 
where the heat generated within the laser rod, due to the 
optical pumping, are removed by a coolant flowing along the 
cylindrical rod surface. As a result of the assumption of 
uniform internal heat generation and uniform cooling along 
the cylindrical surface of an infinitely long rod, the tem- 
perature profiles will be strictly radial. The end effects 
and small variation of coolant temperature in the axial di- 
rection are neglected here. The radial temperature distri- 
bution in a very long cylindrical rod is obtained from the 
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one 


dimensional heat conduction equation: 
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d 2 T , dT 
dr 2 rdr 



[3-D 


where K is the material thermal conductivity and Q is the 
rate at which heat is generated per unit volume. The solu 
tion to equation (3-1) gives the temperature at any point 
along a radius of length r. With the boundary condition 
T(r 0 ) for r=r D , where T(r 0 ) is the temperature at the rod 
surface and r 0 is the radius of the rod, then. 


T(r) = T(r„) * |,(r 2 - r 2 ) 


[3-2] 


This temperature profile is parabolic, with the 
highest temperature being at the center of the rod. The 
heat generated per unit volume, for uniform heat absorb- 
tion, can be expressed as: 


Q - 



(3-3) 


where P a is the total power absorbed by the rod; r Q is the 
rod radius; and L is the length of the rod. The tempera- 
ture difference between the rod surface and the center is 
given by: 

T( 0 ) - Tig = ^ (3-‘ 

The transfer of heat between the rod and the coolant 
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creates a temperature difference between the rod surface 
and the coolant. A steady-state equilibrium temperature 
will be reached when the rate of the total internal genera- 
tion (P a ) of heat is equal to the rate at which heat is re- 
moved from the rod by the coolant. For a typical Nd : YAG 
rod, the equilibrium temperature is reached about 3 1/2 
seconds after laser turn-on. Figure 3-1 depicts an example 
of the radial temperature profile of a 20 pulse per second, 
side-pumped, Nd:YAG laser at thermal equilibrium. 



Figure 3-1. Radial Temperature Distribution 
for a Nd : YAG Crystal 


The heat is generated at the rate of 0.3 Joules per pulse 
and negligible cooling was assumed at the ends of the rod. 
With these conditions, axial thermal gradients do not exist 
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in this rod. 


These temperature gradients generate mechanical radi- 
al, tangential, and axial stresses in the laser rod, since 
the hotter inside area of the rod is constrained from ex- 
pansion by the cooler rod surface. Normally, these stres- 
ses are calculated when one needs to find the birefringence 
and the stress fracture limit of the crystal. For Nd : YAG 
rods, Koechner ^ • 9 has provided a sample of these calcula- 
tions. For example, if the approximate temperature differ- 
ential between the center and the surface of the rod is 
115°C , then the crystal may fracture. The equations for 
these calculations are derived by Timoshenko and 
Goodier . ^ 

Having discussed the basic theory of one dimensional 
temperature distribution and the mechanical stresses caused 
by the thermal gradients, our attention can now be turned 
to the photoelastic effects (birefringence). The mechani- 
cal stresses generate thermal strains in the rod, which in 
turn produce refractive index variations via the photo- 
elastic effect. The refractive index of a medium is speci- 
fied by the indicatrix. In its most general case, the 
indicatrix is in the shape of an ellipsoid. A change of 
refractive index due to strain is given by a small change 
in the shape, size, and orientation of the indicatrix . ^ 
Since Nd : YAG is a cubic crystal, the indicatrix is a 
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sphere. Under stress, the indicatrix becomes an ellipsoid. 
Nd:YAG rods are grown with the cylindrical axes along the 
direction. An illustration of this is depicted in 
Figure 3-2. The light propagates in this direction, and 
thus the change of the refractive index along the [111] 
axis is of interest. 

Since the transverse stresses are in the radial and 
tangential direction (see Figure 3-2), the local indicatrix 
also orients its axis in those directions. In a cylindri- 
cal coordinate system, the photoelastic changes in the re- 
fractive index for the radial and tangential polarizations 
are only obtained after a considerable amount of tensor 
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analysis. This is normally done using the vector Kirchhoff 
diffraction formulation approach. Various treatments of 
these involved birefringence calculations for Nd:YAG are 
covered in detail in references 4, 8 , 12, and 17. The ex- 
pressions for the radial and tangential refractive index 
changes for Nd:YAG are given by®'^ as: 

An, - f C r r* (3-5) 

An, = fV t 3 ' 6 ' 

where 4 n r and fln^ are the radial and tangential stress de- 
pendent changes of the refractive index, respectively; n Q 
is the refractive index of the rod at its ambient tempera- 
ture; a is the thermal coefficient of expansion for 
Nd : YAG ; Q is the heat generated per unit volume; and K is 
the thermal conductivity of Nd:YAG. The coefficients C r 
and C 0 are functions of the el ectrooptical coefficients 
that were determined from tensor analysis for Nd:YAG. Both 
coefficients were calculated by Koechner^ for Nd:YAG. 

Their values are listed as C r = 0.017 and C 0 = -0.0025 re- 
spective! y . 

Having explored the temperature and stresses in a 
Nd : YAG rod, the optical distortions which are a result of 
both the thermal gradients and the stresses can now be dis- 
cussed. The change in the refractive index can be sepa- 
rated into a temperature and a stress-dependent contribu- 
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tion.^ Therefore: 

n[r] = n 0 + An[r] T + An[r] s [3-7] 

where n(r) is the radial variation of the refractive index; 
n Q is the refractive index of the material at ambient temp- 
erature; ^n(r)<ji and 4n(r) s are the temperature and stress- 
dependent changes of the refractive index, respectively. 

The temperature dependent change of the refractive 
index can be expressed as : 2 

An(r] T = [T(r] - T(0]]|jt [3-8] 

where T(r) and T(0) are the temperature at the radius r and 
center of the rod, respectively. With the aid of equa- 
tions ( 3-2 through 3 - 4 ) and equation ( 3 - 8 ), we obtain: * 

An[r] T = - ||r 2 [3-9] 

Comparing equations ( 3 - 5 , 3 - 6 , and 3 - 9 ), we see that the 
refractive index in a laser rod shows a quadratic variation 
with radius. An optical beam propagating along the rod ax- 
is also develops a quadratic spatial phase variation. This 
perturbation is equivalent to the effect of a spherical 
lens-^-4. The focal length of a lens-like medium whose index 
of refraction is assumed to vary according to: 2 

n[r] = n 0 ^ 1 - ^ J [3-10] 
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where b measures the degree 
refraction is given by:^ 

f = 

The total variation of the refractive index is ob- 
tained by introducing equations (3-5, 3-6, and 3-9) into 
(3-7). This gives us : 2 

"(f) = n 0 { 1 - jL § * nJaC,, ]r 2 } (3-12) 

As was previously discussed, the change in the refrac- 
tive index, due to thermal strain, is dependent on the ra- 
dial (r) and tangential (0) stresses where C r = 0.017 and 
C 0 = -0.0025. Therefore, the expression for n(r) in (3-12) 
yeilds two values. Comparing equations (3-10) and (3-12) 
yields the focal length (f') equation : 2 

f’ — — f — + ryT 'l [3 — 131 

1 QU2dT aL r,<t n o J LJ 

The first term in the parentheses is the temperature depen- 
dent variation of the refractive index. The second term 
contains the stress dependent variations of the refractive 
index. Because of the birefringence, NdtYAG has two effec- 
tive focal lengths. One is due to the radial electrooptic 


of variation of the index of 


4n 0 L 


(3-11) 
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(C r ) coefficient and the second is due to the tangential 
(C^) coefficient. 

In our final expression for the focal length of a 
NdrYAG rod, we must also include the effect of rod end dis- 
tortion (bulge). The stresses, that cause bulging, are 
found to occur within a region of approximately one radius 
from the ends for Nd:YAG. 18 This deviation from flatness 
of the rod ends is obtained from: 2 

L(r) = 0 iL o [T(r] - T(0)] (3-H) 

where a is the thermal coefficient of expansion for Nd : YAG 
and L 0 is the length of the end section of the rod where 
expansion occurs. With L 0 = r Q and equation (3-2), we 
obtain : 2 

L(r) = -ar 0 ^- [3-15] 

The contribution to the focal length of the rod caused 
by end-face curvature is given by the thick-lens formula of 
geometric optics-*- 8 as: 

f " ■ 3 e 4i (3 ' 16) 

where the radius of end-face curvature is defined as R = - 
(d 2 L/dr 2 ) -1 . From the expressions (3-14, 3-15, and 3-16) 
follows the focal (f") length of the rod caused by bulge: 2 
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(3-17) 


r = K[aQr 0 (n 0 - 1)] 

Inserting equation (3-17) into (3-13) gives us the final 
expression for the focal length of a Nd:YAG rod:^ 


r _ KAf Jn_ 
1 " P a L 2dT 


aC_ n;? 



(3-18) 


where A is the rod cross-sectional area and P a is the power 
absorbed . 

This expression accounts for the temperature, stress- 
dependent variation of the refractive index, and end-face 
curvature. If one introduces the appropriate material pa- 
rameters for Nd:YAG into equation (3-18), then one finds 
that the temperature dependent variation of the refractive 
index constitutes the major contribution to thermal len- 
sing.^ The stress dependent variation of the refractive 
index and end-face curvature modifies the focal length by 
approximately 20 percent and less than 6 percent, respec- 
tivel y . ^ 

Figure 3-3 depicts the graph of the approximation of 
the focal lengths, using equation (3-18), for a side-pumped 
Nd:YAG rod. The focal lengths were calculated using the 
following parameters: length and radius of the rod was 

7.6cm and 25mm, respectively; dn/dt, the material variation 
of the refractive index was 8.9 x 10“^; index of refraction 
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Figure 3-3. Predicted Thermal Focusing Effects 
in Nd : YAG 

(n 0 ) was 1.825; thermal conductivity (K) was 0.103; a, the 
thermal coefficient for expansion was 7.9 x 10“^; the radi- 
al (C r ) and tangential (C^) variation of the refractive in- 
dex were 0.017 and -0.0025, respectively; and, the power 
absorbed (P a ) by the rod (as heat), was 6.119w/cm^. 

This completes the one dimensional theory of thermal 
lensing for a Nd:YAG rod. The same analysis holds for 
Ti:AL203 and other laser materials. As was previously 
stated, the temperature dependent variation of the refrac- 
tive index will be the only aspect implemented in the Ther- 
mal Lensing Model. The implementation of birefringence and 
bulge effects will be left as future enhancements to the 
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model. In the next section, the formulation of the Thermal 
Lensing Model is developed. 

C. Two Dimensional Model Formulation. 

Formulation of the thermal lensing ray-trace model 
requires basic understanding of various optical theories. 
These include the theories of ray tracing 16 , wave optics in 
crystals 11 ' 19 , Fresnel diffraction theory of waves in lens- 
like mediums 1 ^, geometrical -optics 16 , etc. It was not the 
purpose of this thesis to become trapped into the finer de- 
tails of the referenced optic theories. Nor do we intend 
to repeat the derivations of the referenced theories. If 
more detailed information is desired, then we refer the 
reader to the attached list of references and other litera- 
ture on the subject. 

The basic task of this research is to formulate a com- 
puter ray-trace model that allows for a first approximation 
of focal length. The model is to be based on the tempera- 
ture variation of the refractive index. The approximations 
of focal length are then obtained on specific crystals 
operated with specified parameters. This will be done by 
modifying the focal length expression (3-11) for a thin 
lens. The focal lengths are then approximated, by calcu- 
lating the phase shift in the index of refraction, as the 
different rays of a plane wave front are tracked through a 
lens-like crystal medium. The localized change in the ray 
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direction as it transverses a lens-like medium is assumed 
to obey Snell's law which is given by: 19 

rifSinOp = i^sinOj [3-19] 


where n^ and n 2 are the index of refractions of two adja- 
cent mediums, respectively; and 0^ and 0 O are the input and 
output angles of the ray at the boundary between the adja- 
cent mediums. 

To assist in the formulation of the model, it is ap- 
propriate to introduce, and explain, two tables and a fig- 
ure before proceeding. Table 3-1 depicts an illustration 
of the output block of radial and axial temperature pro- 
files, provided by the Temperature Distribution Model. 

This output is based on a 20pps, side-pumped, Nd : YAG laser. 
The laser was operated with moderate side cooling and rod 


Table 3-1. Block of Radial and Axial Temperature 
Profiles for a Nd:YAG Crystal 



Temperature Distribution 

Radial Diffusivity 

Radius of Rod (cm) 

Length of Rod (cm) 

Time ( sec ) 

26.95 
0 . 25 
7.6 
3.5 


R\L 

0.00 

0.50 

* . . 

1.00 

0.00 

2 . 9139 

2.9139 


2 . 9139 

0 . 20 

2.8020 

2.8020 


2 . 8020 

0.40 

2 . 4660 

2.4660 


2 . 4660 

0.60 

1.9047 

1.9047 


1.9047 

0 . 80 

1.1276 

1.1276 


1.1276 

1.00 

0.2248 

0.2248 


0 . 2248 
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ends were not. cooled. In the vertical direction, are the 
radial temperatures starting from the center of the rod 
(top) and increments outward to the surface of the rod 
(bottom). In the horizontal direction, are the axial tem- 
peratures starting from the front face (left) of the rod 
and increments to the back end (right) of the rod. In this 
particular example, a set of 10 equally spaced increments 
was chosen to represent the 25mm x 7.6cm rod in both the 
axial and radial directions. Additionally, in Table 3 1, 
the temperatures are observed only to vary in the radial 
direction. The displayed block of temperatures in Table 
3-1 is an example of a side-pumped Nd:YAG rod at thermal 
equilibrium, when rod ends aren't cooled. 

Table 3-2 displays a block of radial and axial index 
of refractions for the temperatures presented in Table 3-1. 


Table 3-2. Block of Radial and Axial Index of 
Refraction for the Temperatures in 
Table 3-1 


Index of Refractions 


R\L 


0 .00 


0.00 
0.20 
0.40 
0 . 60 
0.80 
1.00 


1.825026 

1.825025 

1.825022 

1.825017 

1.825010 

1.825002 


0.50 


1.825026 
1.825025 
1 . 825022 
1.825017 
1.825010 
1.825002 


1.00 


1.825026 

1.825025 

1.825022 

1.825017 

1.825010 

1.825002 
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The block of temperature variations of the Index of Refrac- 
tions were calculated from the expression: 

n(r,z) = n 0 + aT(r,z] (3-20] 

where n Q is the index of refraction for the rod at its am- 
bient temperature and a is the dn/dt specified for the rod. 

Table 3-2, shows that the index of refraction is only- 
changing in the radial direction which was expected for a 
side-pumped Nd:YAG rod with no end cooling employed. For 
end-pumped rods, the index of refraction would change in 
both the radial and axial direction. The incremental par- 
tition of the index of refractions in Table 3-2 are the 
same as the partitions of temperatures displayed in Table 
3-1. 

The user specified partition, selected for use in the 
Temperature Distribution Model , governs the partition of 
the Thermal Lensing Model. The partition of 10 equally 
spaced increments are normally chosen because it is more 
than adequate for first approximations. Increasing the 
number of partitions does increase the accuracy of the ap- 
proximation. However, the time required to run the Tem- 
perature Distribution Model increases linearly as the num- 
ber of partitions increases. Comparison of the temperature 
outputs for a partition of 10 and 20 equally spaced incre- 
ments showed a mere 0.4 percent difference in the approxi- 
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mated temperatures. 

Figure 3-4 illustrates the three types of pumping beam 
profiles used to pump a laser rod. All three types of pum- 
ping profiles can be used to end pump a rod. (Normally, 
only the circular and parabolic beam profiles are used to 
describe the distribution of heat deposition in side-pumped 
operations . ) 



Figure 3-4. Types of Pump Beam Profiles 


The above mentioned types of pump beam profiles have 
been implemented in the Temperature Distribution Model. In 
the thermal lensing ray-trace model, the focal length of 
the rod irrespective of the type of pumping beam used can 
be approximated. However, only the circular pump beam has 
been implemented for distortion calculations. Implementa- 
tion of distortion calculations for parabolic and Gaussian 
beam profiles have been left as a future model enhancement. 
Table 3-1 represents a matrix of temperatures for a specif- 


37 



ic rod operating with user specified parameters. It serves 
as one of the primary inputs to the Thermal Lensing Model . 

A radially dependent phase delay is experienced by a 
wave front propagating through a medium with a varying in- 
dex of refraction n. Thus, a thin slab of the crystal me- 
dium can be approximated as a thin lens. The lensing model 
traces wave fronts propagating through the rod starting 
from several radial positions at the rod front face. As 
the ray propagates through the rod, a series of thin lenses 
are setup and the temperature dependence of the index of 
refraction is calculated at specified locations. Using 
these calculations, a focal length associated with each 
tracing ray can be estimated. An average focal length for 
the rod is obtained by taking the mean of the focal lengths 
of the individual tracing rays. The specified axial incre- 
ment of the temperature matrix determines the number of eq- 
uivalent thin lens placed in series that the tracing ray 
must pass through. For example, with an axial incremental 
spacing of 0.1, as depicted in Table 3-1, 20 equivalent 
thin lens are placed in series from the front to the back 
of the rod. 

A plane wave propagating through a lens-like medium 
acts as a spherical thin lens and diverges. Using this 
fact, the thin lens focal expression (3-11) was changed to: 


trr 0 ) 2 

AzL ]>X(r,z) - njO.z)] 


(3-21) 
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where r Q and r are the exiting relative radius of the ray 
and the radius of the rod in cm, respectively; z is the 
specified axial increment of the rod (dimensionless); L is 
the length of the rod in cm; and n r (r,z), n o (0,z) are the 
temperature dependent index of refraction for the radial 
point of interest and for the center of the rod, respec- 
tively. The focal length expression (3—21) can be applied 
upon determination of the quantities r Q/ n r (r,z) , and 
n o (0,z). Calculating these scalar quantities are the es- 
sence of the thermal lensing ray-trace model. 

The radial variation in temperature, as confirmed by 
use of curve fitting algorithms, is quadratic in nature. 
This fact allowed the selection and development of a bi- 
nomial interpolation algorithm to obtain values of tem- 
perature at points other than the selected matrix points. 
The binomial interpolation algorithm is based on the ap- 
proximation : 

ax^ + bx + c = T (3-22) 

where a,b,and c are constant coefficients, x is either the 
radial (r) or axial (z) coordinate, and T is the tempera- 
ture. By invoking this algorithm using a set of three col- 
linear temperature (input) values along with their coordin- 
ate specifications, we can obtain a value for the tempera- 
ture at any other point not specified in the temperature 
matrix . 
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Figure 3-5. Partial Illustration of the Thermal 
Lens i n g Matrix With a Series of Thin 
Lens Overlaid 


The temperature matrix (obtained from the input file) 
is illustrated in Figure 3-5. A wave front ray will be 
traced starting at the r^ + p radial input position through 
a series of equally spaced ( Az/2 ) thin lenses. The r j_ + ]_ ray 
was selected to illustrate the definition of the input and 
output angles that occur as the ray passes through a thin 
lens. These angles are necessary to calculate the exiting 
radius of the ray (r 0 ) for use in the focal length expres- 
sion (3-21) and distortion calculations. 

The ry ray (center of the rod) is also traced to ob- 
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tain a summation of the temperature variation of the index 
of refraction at the center of the rod, n o (0,z). As illus- 
trated in Figure 3-5, temperatures through T 9 are ob- 
tained from the temperature matrix provided as input. In- 
voking the binomial interpolation algorithm allowed the 
calculation of the intermediate temperatures T^ c through 
Tg c for use in the thermal lensing scheme. This completes 
the basic setup of the thermal lensing scheme and provides 
the framework to proceed with the ray-tracing of specific 
input rays . 

A collimated plane wave entering the rod will have an 
input angle ( 8 ^ ) equal to zero. As it passes through the 
first lens, the angle (9 m ) between the normal to the sur- 
face of constant temperature and the rod axis will be cal- 
culated at A z/2. To calculate 8 ^, we determine first the 
surface of constant temperature by obtaining the r and z 
derivatives of expression (3-22). The derivatives are giv- 
en by : 

= 2av i + b (3-23) 

= 2a[^j + b (3-24) 

where r^ is the input radius, A z is the axial incremental 
spacing, and a,b are the calculated coefficients using the 
appropriate radial and axial temperature points in the in- 
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terpolation algorithm. The angle 9^ at Azj7 is then, calcu- 
lated using the trigonometric expression: 

9 » = ,an ‘(f ] (3 - 25) 

On finding the ang 1 e ' the next step is to determine the 
specific values of the index of refraction on each side of 
the imagined lens. Due to the fact that the temperature 
dependent index of refraction changes very little between 
adjacent points (<< 1), it was decided to compute the index 
of refraction on each side of 9^ according to the following 
expressions : 

n, = n 0 + — * 2 — J [3-26] 

n 2 = n 0 + J P-2 7] 

where n Q is the index of refraction of the rod at its am- 
bient temperature; a is the specified dn/dT of the rod; and 

T 2 ' t 5 / an< * Tj 4 (same as T 2 C in the first stage of tracing) 

are the specific temperature points before, at, and after 
the position of the imagined lens. 

By invoking Snell's law, the output angle (0 O ) of the 
ray passing through the thin lens can be obtained by sub- 
stitution in the expression: 
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(3-28) 


0 O = 0 m - sin_1 [57 sin ( 0 m - 0 ij] 

After finding the output angle, the change in output radius 
( r 0 ) can be calculated from the trigonometric expression: 

r 0 = ytan 0 O + r m (3-29) 


where r m is the midpoint radial position of the ray (same 
as r-L + i for initial entry of the ray into the rod). The 
output angle (0 O ) and output radius (r Q ) are then sent to a 
recursive numerical algorithm to refine the output values 
at A z to meet operator specified accuracy criteria for 0 O . 

The phase shift in the temperature variation of the 
index of refraction can be calculated, for the initial 
entry of the ray into the rod, and is proportional to the 
expression : 


n r (r,z) 



(3-30) 


and for later increments of the rod to: 


rV(r.z) = n, 


"rjlr 


1 ♦ ^ + n 2 1 ♦ ^ 

N (W N (W 


r J r 


(3-31) 


43 


The individual calculations from each lens, as the ray pro- 
pagates through the rod, are retained and summed to deter- 
mine the value for n r ( r , z ) which is then used in the focal 
length expression (3-21). 

One-to one distortion mapping of a ray propagating 
through a lens-like medium is a straight forward calcula- 
tion when dealing with side and end-pumped circular beams. 
The input irradiance (I) of the collimated plane wave is 
constant as it enters the crystal medium. The incremental 
change in irradiance is then calculated at the output end 
of the rod from the expression: 


. - r . 2 W 

r|(i+l] - r|Ci) 


[3-32] 


where r Q , r ^ are the output and input radius of the ray be- 
ing traced, respectively. The percentage of total distor- 
tion (D), for each ray, is then calculated from the ex- 
pression : 


D = 100 


f r|[i+l) - iffi+1) ''l 

L J 


[3-33] 


This completes the formulation of the thermal lensing 
ray-trace model. Now, it’s just a formality of computer 
record keeping to invoke the binomial interpolation 
algorithm, using r Q , to calculate T^q, T^, and T 12 , reset 
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the output radius ( r c ) of the previous rod segment equal to 
the input radius ( r^ ) of the next rod segment, reset the 
output angle (0 O ) equal to the input angle (8^), and incre- 
ment forward to the next imagined thin lens. The calcula- 
tions are repeated until the ray exits the rod. When the 
ray exits the rod, the focal length and distortion associ- 
ated with the ray are calculated and retained. Then, the 
ray-trace program increments upward in the matrix of tem- 
peratures to trace the next input ray. Upon completion of 
the ray tracing for a block of temperatures, the mean of 
the individual focal lengths and distortion present for 
that block of temperatures are then calculated. The model 
then increments forward to the next block of input tempera- 
tures (at a later time) and repeats the process until the 
input file of temperatures are exhausted and a time depen- 
dent average focal length are obtained. 

The Thermal Lensing Model program code is contained in 
Appendix A. An example of the lensing model inputs and 
output are enclosed in Appendix B and C, respectively. 
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CHAPTER 4 


Analysis and Discussion of Thermal Lensing 

The next generation of solid-state lasers (used for 
remote sensing of the atmosphere from space platforms), are 
under development. Various new lasing materials and laser 
conf iguations are being investigated. Currently, the so- 
lid-state, end-pumped diode lasers are the most promising 
candidates due to their small size, reliability, perform- 
ance, low power and cooling requirements, etc. Analyzing 
thermal lensing effects for these end-pumped conf iguations , 
using various cooling schemes are of interest to the de- 
signers of laser systems. Theoretical formulations, numer- 
ical calculations and experimental data do not exist for 
end-pumped conf iguations . 

Before proceeding with the analysis of thermal lensing 
in end-pumped Nd:YAG and Ti:AL203 laser rods, an assessment 
of the validity of the Thermal Lensing Model was conducted. 
The validation was based on the comparison of model results 
with the experimental data for a repetitively pulsed side- 
pumped Nd:YAG laser. This case was interesting from two 
points. It was the only laser found, with actual experi- 
mental data on thermal focal length. The second point of 
interest is that the laser was side-pumped. The later 
point allowed for the extrapolation of Koechner's predic- 
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tion which are based on a CW side-pumped Nd : YAG laser. 

Thus, the model results were compared to both the experi- 
mental results and Koechner's theoretical formulation for a 
long continuously side-pumped rod. 

A. Model Validation. 

Research of the literature shows that theoretical for- 
mulation and experimental data on thermal lensing effects 
do not exist for repetitive pulsed end-pumped lasers due to 
their recent development. In fact, very little data exist 
on repetitively pulsed side-pumped lasers. The vast major- 
ity of data published such as theoretical formulations, 
theoretical calculations, experimental results, etc. are 
based on CW side-pumped lasers dating back to the early 
1970's. Most of these early analyses were conducted by Ko- 
echner 5 ' 9 and are still being referenced to as the authori- 
tative source on the subject today. 

A recent experiment, depicting the results of thermal 
focusing for a repetitive pulsed side-pumped Nd : YAG laser 
has been completed. ^0 The laser rod used in the experiment 
had a length of 7.6cm, a radius of 25mm, was side-pumped 
uniformly, and cooled along its length (ends were not 
cooled). Each pump pulse had a duration of 130us and an 
energy level of 10 Joules. Experimental calculations indi 
cated that approximately 3.5 percent of the pump energy was 
absorbed into the rod as heat. 
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Both, the Temperature Distribution Model and Thermal 
Lensing Model were executed using the Nd:YAG material char- 
acteristics and the laser configuation of the experimental 
setup. The modeling and experimental results, along with 
Koechner's prediction of focal length from equation (3-18) 
are shown in Figure 4-1. (It should be noted that Koech- 
ner's prediction, used in Figure 4-1, did not include bire- 
fringence and bulge effects.) 



Figure 4-1. Thermal Focus ing in Nd:YAG Rod 


Figure 4-1 shows that the model's prediction of the 
thermal focal length are in good agreement with the experi- 
mental results and Koechner's theoretical prediction of fo- 
cal length for an equivalent power problem. Thus, the ray- 
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trace model, based on the thermal gradients provided by the 
Temperature Distribution Model, will give a fairly accurate 
first order approximation. 

Transient thermal focal length calculations of a side- 
pumped Nd:YAG rod do indicate an appreciable deviation from 
the standard steady-state power calculations using Koech- 
ner's prediction equation (3-18). Starting from an inac- 
tive condition, the Nd : YAG rod thermal focal length did 
stabilize at about 3.5 seconds. However, the steady-state 
calculations from the model at the end of each pump pulse 
indicates a higher focusing power of about 15-25 percent at 
frequencies of 4-20 hertz over the results as predicted by 
Koechner. This was not unexpected as the model was exe- 
cuted in an repetitive pulsed mode. Where as, Koechner' s 
prediction are based on a thermally stabilized CW pulsed 
(same power) rod. 

The transient side-pumped thermal focal length varies 
from 28 meters at 100ms to 6.97 meters at the steady-state 
equilibrium condition. The time dependent focal length, as 
shown in Figure 4-2, was computed at the end of each pump 
pulse for the above described laser rod. The rod was pum- 
ped by a 20 Hz, 10 Joule per pulse flash lamp. The hori- 
zontal line was obtained by using Koechner's formula for a 
thermally stabilized CW pulsed (same power) rod whose sur- 
face temperature was kept at a fixed value. 
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Figure 4-2. Transient Thermal Focal Length for 
a Side-Pumped Nd : YAG Rod 


Thermally induced focusing varies significantly in the 
time interval between pulses. Even at a 10Hz repetition 
frequency, the variation in the Nd : YAG thermally induced 
focal length was approximately 10 percent. 

The model was also used to estimate the impact of sur- 
face cooling on the focal length. Focal length was found 
to be very sensitive in the low end of the surface cooling 
parameter. Figure 4-3 displays the results of focal length 
versus cooling in side-pumped Nd : YAG rods. The cooling 
parameter is a dimensionless quantity obtained by multi- 
plying the length of the rod by the surface heat transfer 


50 




1 /FOCAL LENGTH (DIOPTERS) 


coefficient and dividing by the bulk thermal conductivity 
of the rod material. Koechner's formulation, by compari- 
son, is for a fixed surface temperature. 
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Figure 4-3. Focal Length vs Cooling i n Side- 
Pumped Nd:YAG Rod 
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Figure 4-3 is also important to system designs. It 
shows that above a certain point, very little is gained 
from additional cooling. Therefore, the most efficient 
cooling for the laser is critical for spaced based applica- 
tions. Power consumption in space platforms are one of the 
most critical design parameters. Additional unneeded 
cooling may strain the limited power resources available in 
space . 
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One-to-one mapping of intensity in the repetitive pul- 
sed side-pumped Nd:YAG laser showed distortion to be negli- 
gible during the rod transient period. This result was 
expected due to the presence of only radial thermal gradi- 
ents (no axial gradients) as a result of the uniform pum- 
ping of the rod and the absence of end-cooling. 

B. Analysis of Thermal Lensing in End-Pumped Nd : YAG and 

Ti:AL203 Laser Rods. 

Thermal lensing calculations have been obtained for 
the two end-pumped laser rod conf iguations depicted in Ta- 
ble 4-1 below. Both of these laser crystals are currently 
undergoing feasibility experiments in the laboratory at 
NASA Langley Research Center. It should be noted here that 
the comparison is not between crystals of different dimen- 
sions, but between cyystals of practical dimensions. Both 

Table 4-1. End-Pumped Laser Rod Characteristics 


Parameters 

Nd : YAG 

Ti : AL2O3 

Rod Length (cm) 

1.0 

0 . 6 

Rod Radius (cm) 

0 . 4 

0 . 2 

Thermal Diffusivity (cm 2 /s) 

26.947 

8.87 

Thermal Conductivity (W/cm °K) 

0.103 

0.330 

Material dn/dt 

8. 9xl0~ 6 

1 . 28xl0~ 5 

Index of Refraction 

1.825 

1.7654 

Optical Attenuation 

3.0 

3.0 

Pump Wavelength (microns) 

1.15 

0 . 6438 

Pump Energy Absorbed (J/cm^) 

0.050 

0.050 

Pump Pulse Duration (seconds) 

0.010 

0.010 


crystals were end-pumped by a Gaussian beam where the ratio 
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of the beamwidth to the rod radius is 0.5. With the excep- 
tion of the rod material characteristics, all other parame- 
ters such as side and end cooling, energy absorbed, pump 
pulse duration, etc. were the same for both rods. Both, 
the Temperature Distribution Model and the Thermal Lensing 
Model were executed using the appropriate material charac- 
teristics of the rod and the laser configuation depicted in 
Tabl e 4-1 . 

Figures 4-4 and 4-5 presents the single pulse tran- 
sient behavior of thermal focal length for Nd:YAG and 
T.i:AL 2 C >3 crystals, respectively. Both rods were end-pumped 



o.oi 

TIME (SEC) 


Figure 4-4. Sin g le Pulse Transient of Thermal 
Focal Length i n Nd:YAG Cry s tal 


with a single pulse of 10ms duration. The energy absorbed 
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Figure 4-5. Single Pulse Tra n sient of Thermal 
Focal Length in Ti : ALoOA._Cr yst . a l 


into the rod as heat was assumed to be 50mJ per pulse. Mo- 
derate side cooling parameter of 20 was employed. Ends 
were not cooled. (The cooling parameter is a dimensionless 
quantity obtained by multiplying the length cf the rod by 
the surface heat transfer coefficient and dividing by the 
bulk thermal conductivity of the rod material.) Starting 
from an inactive condition, the thermal focal length for 
Nd : YAG ranged from 205 meters at 2ms to the minimum of 44.8 
meters at 10ms when the pump pulse was turned off. For 
Ti : AL 2 0 3 , the thermal focal length ranged from 18 meters at 
1ms to 2.87 meters at 10ms. Although some residual heat 
still exist in the rods, as the rods relaxes to its ambient 
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condition, the thermal focal lengths are no longer measur- 
able 1.30 seconds and 250ms after pump shutdown for Nd : YAG 
and Ti : AL 2 O 3 rods, respectively. For all practical pur- 
poses, as far as focal length is concerned, the rods have 
returned to their ambient state at this point. 

The transient thermal focal length profiles for Nd:YAG 
and Ti:AL 2 03 during multiple pulse operations are shown in 
Figures 4-6 and 4-7. As seen in Figure 4-6, the Nd : YAG rod 



Figure 4-6. Transient Thermal Focal Length for 

a Multiple Pulsed, End-Pump Nd:YAG Rod 


thermal focal length stabilized at about 3.5 seconds. For 
Ti:AL 203 , the rod thermal focal length stabilized at about 
300ms. In both rods, it was found that the thermal focal 
lengths stabilized at the same time the temperatures became 
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Figure 4-7. Transient Thermal Focal length for a 

Multiple Pulsed, End-pumped TirALoO j Rod 


stabilized which was expected. The Nd : YAG transient end- 
pumped thermal focal length varied from 44 meters at 10ms 
to 6.11 meters at the steady-state equilibrium condition. 
For Ti:AL 203 , the thermal focal length varied from 2.87 me- 
ters at 10ms to 2.66 meters at steady-state equilibrium. 
These time dependent focal lengths, shown in the solid 
lines, in Figures 4-6 and 4-7 were computed at the end of 
each pump pulse for the described laser rods. Both rods 
were pumped at a 20 Hz rate. Energy absorbed into the rods 
as heat was assumed to be 50 mJ per pulse. 

Thermally induced focusing varied significantly in the 
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time interval between pulses as shown by the dashed lines 
in Figures 4-6 and 4-7. The variation was approximately 26 
and 75 percent for Nd : YAG and TiiAI^C^, respectively. The 
time dependent focal length between pulses were computed at 
the start of each pump pulse for the above rods. Not ac- 
counting for the temporal variation of the thermally in- 
duced focal length during the pulse relaxation period, may 
lead to significant errors during design. 

The next factor looked at was thermal focusing versus 
pump frequency and the impact of side and end cooling on 
thermal focal lengths. Figures 4-8 and 4-9 depict the re- 
sults for Nd : YAG and Ti:AL 203 , respectively. Pump frequen- 



Figure 4-8. Thermal Focusing in End- 
Pumped Nd : YAG Rod 
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cy was varied between 1 and 20 HZ. For each frequency, the 
rod was allowed to stabilize (3.5s for Nd:YAG and 300ms for 
TirAL^C^) before calculations were made. In both figures, 
moderate side cooling parameter of 20 was employed. For 
calculations involving both side and end cooling, end 
cooling was set equal to the side cooling. As shown in 
Figure 4-8, the impact of end cooling caused a significant 
decrease in the Nd:YAG thermal focal length above 2 Hz. 

This result was expected as the radial and axial tempera- 
ture gradients near the rod ends become smoother due to end 
cool ing r 
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For Ti : AL 2 0 3 , very little impact from end cooling was 
experienced until the 12 HZ pump frequency was reached. 
However, in comparison to Nd : YAG where the thermal focal 
length was significantly decreased, the opposite effect was 
obtained for Ti:AL 2 03 . End cooling of Ti:AL 2 0 3 crystals 
causes small increases in the thermal focal length as the 
pump frequency is increased. This result was not expected. 
However, upon further analysis of the temperature profiles 
for Ti : AL 2 0 3 , it was apparent that the radial and axial 
temperature gradients near the ends were greater due to end 
cool ing . 

Finally, the model was used to estimate the impact of 
surface cooling on the induced thermal focal lengths for 
end-pumped Nd : YAG and Ti:AL 2 0 3 laser crystals. The rods 
were pumped at a 20 HZ rate and energy absorbed into the 
rod as heat was assumed at 50mJ per pulse. Figures 4-10 
and 4-11 presents the results of focal length versus side 
cooling for these conf iguations . The focal length was 
found to be very sensitive only on the low end of the sur- 
face cooling parameter for both rods. Therefore, sizing 
end-pumped lasers for correct cooling is critical for sp- 
aced based applications as previously explained. 

Distortion calculations were not performed on these 
end-pumped lasers as the function for mapping of Gaussian 
beams has not been implemented in the model (left as a fu- 
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Figure 4-10. Focal Length vs Side Cooling 

in End-Pumped Nd:YAG Rod 



Figure 4-11. Focal Length vs Side Cooling in 
End-Pumped TirALoC^ Rod 
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ture model enhancement). Unlike the side-pumped configu- 
ation, previously discussed, where distortion was non- 


existent at thermal 
are that distortion 
on the order of 3-5 


equilibrium, preliminary indications 
in the end-pumped conf iguations will be 
percent at thermal equilibrium. 
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CHAPTER 5 


Summary and Conclusions 

An analytic expression for approximating the time-de- 
pendent thermal focal length in finite solid-state laser 
rods has been derived. The analysis is based on the tem- 
perature variation of the material refractive index caused 
by optical pumping of these rods. This analysis is also 
applicable to situations where the laser rod is side or 
end-pumped by circular, parabolic, or Guassian beam pro- 
files. Short or long pulse, continuous wave (CW), or re- 
petitive pulsed operations are also applicable. The repet- 
itive pulsed, end-pumped solid-state diode laser configu- 
ations are of current research interest as the next gener- 
ation of lasers will come from this category. 

Several quantities were found to be relevant to this 
analysis. These quantities were the specific thermal pro- 
files of the rods, type of optical pumping employed, type 
of cooling scheme employed (side and end-cooling parame- 
ters), and the specific material characteristics of the 
rods . 

The Thermal Lensing Model was formulated using geomet- 
ric ray tracing based on the focal length expression for a 
thin lens. The focal lengths are then approximated, by 
calculating the phase shift in the index of refraction, as 
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the plane wave is tracked through a lens-like crystal medi- 
um. 

Thermally induced temperature profiles, for end-pumped 
laser rods operating with specified parameters, are approx- 
imated by the Temperature Distribution Model . These ther- 
mal profiles are then inputted to the developed Thermal 
Lensing Model . The Thermal Lensing Model is then used to 
approximate the transient thermal focal length and estimate 
the final rise in thermal focal length for the laser rod 
under evaluation. Various pumping and cooling schemes can 
be evaluated for their effect on the induced thermal focal 
length. For circular pump beam profiles, distortion calcu- 
lations are also approximated by the model . 

The numerical examples for end-pumped, repetitively 
pulsed Nd : YAG and Ti:AL 203 laser rods showed the following 
major results. 

a. The minimum induced thermal focal length is 
achieved at the time of pump pulse turn off 
for both rods operating under single long 
pulse conditions. 

b. For multiple pulse operations, the minimum 
induced thermal focal length was achieved when 
the rod reached its thermal equilibrium cond- 
ition (3.5 seconds for Nd:YAG and 300ms for 

Ti : AL 2 0 3 ) . 
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c. Thermally induced focusing varied significant- 
ly in the time interval between pulses (26 
percent for Nd:YAG and 75 percent for Ti : - 
AL 2 0 3 ) . 

d. Employment of end-cooling significantly de- 
creased thermal focusing in Nd:YAG above the 2 
Hz pump frequency. For Ti:AL 2 C> 3 , the opposite 
effect was found. Small increases were 
noticed in the thermal focal length for pump 
frequencies above 12 Hz. 

e. Under conditions when only the side cooling 
was varied, the induced thermal focal length 
was found to be very sensitive at the low end 
of the scale for both rods. 

The model was validated against experimental data and 
theoretical calculations of thermal focal length for a re- 
petitive pulsed, side-pumped NdtYAG laser. It is seen that 
the prediction of thermal focal length is in good quantita- 
tive agreement with experimentally observed behavior. Thus, 
the Thermal Lensing Model, based on the thermal profiles 
provided by the Temperature Distribution Model, will give a 
fairly accurate first order approximation of the thermal 
focal length. 
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Appendix A 


Computer Program Listing 

1. Functional Description 

The Thermal Lensing Model (THERLENS . FOR) is not a 
standalone computer program. Its primary input is depen- 
dent on the output of a previously developed theoretical 
model titled: Temperature Distribution Model for finite op- 
tically pumped laser rods. The output of the Temperature 
Distribution Model (blocks of temperature data versus 
time), along with other thermal lensing constants (see Ap- 
pendix B), are inputted to the lensing model. The model 
then calculates the finite changes in focal length and dis- 
tortion effects due to the temperature gradients within the 
laser rod. Appendix C depicts the Thermal Lensing Model 
output . 

Currently, the lensing model can calculate the thermal 
lensing effects for side and end-pump configurations and 
distortion approximations for circular pump beams only. 
(Note: Distortion calculations for parabolic and Gaussian 
pump beams will be implemented in the future.) The optical 
axis of the laser rod is assumed to be along the axis of 
the rod. 

2. Environmental Characteristics 
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This program was written in the Fortran language. It 
was compiled and linked using the Microsoft Fortran 77 
Optimizing Compiler Version 4.10. The Thermal Lensing Mod- 
el can be executed on any IBM/IBM compatible 286/386AT mi- 
cro-computer with a 287/387 Math Coprocessor. This version 
can also be executed on the VAX 750 mini -computer . The 
source program size on the AT-machine is 46,527 bytes and 
the executable program is 69,621 bytes. The enclosed pro- 
gram code contains all changes through March 1, 1990. 

3. Thermal Lensing Program Listing 

(Program listing starts on next page) 
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Thermal Lensing Program Listing 


03-01-90 

18:42:07 


Source Line 


Microsoft FORTRAN Optimizing Compiler Version 4.10 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 


1 
2 

3 

4 

5 

6 

7 

8 
9 

10 
11 
12 

13 

14 

15 

16 

17 C 

18 C 

19 C 

20 C 

21 C 

22 C 

23 C 

24 C 

25 C 

26 C 

27 C 

28 C 

29 C 

30 C 

31 C 

32 C 

33 C 

34 C 

35 C 

36 C 

37 C 

38 C 

39 C 

40 C 

41 C 

42 C 

43 C 


$DEBUG 
C 


A COMPUTER PROGRAM SIMULATION MODEL TO COMPUTE THERMAL LENSING FROM 
TEMPERATURE DISTRIBUTION PROFILES IN OPTICALLY PUMPED LASER CRYSTAL 
RODS, (YY.FOR) 

******************************************************************* 


THERLENS.FOR: Version 1.0 

Language: Microsoft 


Fortran77 Version 4.10 


Written by: Vincent G. 


Brackett 

Graduate Student Applied Mathematics 
Hampton University 
Hampton, Va. 23668 


Thesis Advisor: 


Dr. Usamah O. Farrukh 

Department of Engineering 
Hampton University 


For: 


GRANT: 

Monitor: 


NASA Langley Research Center 
Hampton, Va. 23666 

NAG 1-666 
Philip Brockman 


Date: Last Revised: March 1, 1990 

******************************************************************* 


THIS PROGRAM USES THE TEMPERATURE OUTPUT PROFILES OF A PREVIOUSLY 
DEVELOPED COMPUTER SIMULATION MODEL (THAT MODEL PREDICTS THE 
NORMALIZED TEMPERATURE DISTRIBUTIONS VERSUS TIME IN FINITE SIDE AND 
END-PUMPED LASER CRYSTAL RODS) TO PREDICT THERMAL LENSING (FINITE 
CHANGE IN FOCAL LENGTHS) CAUSED BY THE RESULTANT TEMP GRADIENTS IN 
THOSE RODS. THERMAL DISTORTION OF THE PUMP BEAM ARE CALC AND AVG 
VALUES OF THE FOCAL LENGTHS AND DISTORTION ARE PRESENTED. 

CAUTION: The matrix mesh of the Temperature Distribution Model and 

this model must be the same size. To change matrix mesh 
size, only the parameter statement needs to be changed. 


67 



44 

45 C 

46 

47 

48 C 

49 

50 C 

51 

52 

53 C 

54 

55 C 

56 

57 

58 

59 

60 C 

61 
62 

63 

64 C 

65 C 

66 C 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 


PARAMETER (N=11,M= 11} 

DIMENSION TEMP(N,M) , INDEX(N,M) ,RAD(M) ,LON(N) ,DISTOR(M) 

DIMENSION ANGI(M) ,ANGM(M) ,ANGO(M) ,REND(M) , FOCAL (M) 

REAL LAMBDA , L AMBDL , K , LRMIN , LRMAX , LR , LENG , LON , IOR , INDEX , N1 , N2 , LONI 

DOUBLE PRECISION PI , SDMIOR , SDMINT , SIORO , SIORJ , SINTJ , ANGI J , ANGOJ 
DOUBLE PRECISION RI,RM,RO,Z,A,B,C 

CHARACTER* 1 XX (140) 

COMMON / BLOCK1 / IOF , DT , ACC , N1 , K2 , ZI , RI , RM,TM, LONI , I 
COMMON /BL0CK2/R1,R2,P.3,Z1,Z2,Z3 
COMMON /BLOCK3/T1 , T2 ,T3,T4 ,T5 ,T6 ,T7 ,T8 ,T9 
COMMON /BLOCK4/T1C ,T2C ) T3C,T4C, T5C , T6C r T10 , Til , T12 

OPEN ( 4 , FILE= ' THERLER1 , INP ' , STATUS: ' OLD ' ) 

OP EN { 5 , FI LE= ' THERL EN2 . 1 NP ' , STATUS? 1 OLD ' ) . 

OPEN( 6 , FILE= ' THERLENS , OUT ' , STATUS^ 'UNKNOWN ' ) 

...{PROCESSING OF INPUT FILE ONE)... 


WRITE(6,*)' ' 

WRITE(6,*)‘ THERMAL LENSING OUTPUT FILE 

WRITE(6,*)' ' 

WRITE(6,*) ' INPUT FILE *1 

WRITE(6,*)' 


READ(4, ' (80A1) ' ) (XX(I) ,1=1,80) 

WRITE{ 6 , * ) <XX(I) ,1=1 ,80) 

READ(4, * { 80 A1 ) ' ) (XX ( I ) , 1 = 1 , 80 ) 

WRITE(6,*) { XX ( I ) ,1=1,80) 

READ(4,*) WL 

WRITE(6,*) '..WAVELENGTH OF INTEREST IN UM..(WL).: ’ ,WL 
READ { 4 , * ) TOI 

WRITE( 6 ,*) ' . .TEMPERATURE OF INTEREST IN DEG C..(TOI).: \TOI 
READ{4,*) IOR 

WRITE ( 6 , *) '..INDEX OF REFRACTION. . (IOR) . : ',IOR 
READ (4,*) DT 

WRITE(6,*) '..ALPHA (dn/dt) . . (DT) . : \DT 
READ(4,*) TNF 

WRITE (6,*) '..TEMP NORMALIZATION FACTOR.. (TNF).: 1 ,TNF 
READ(4,*) TINF 

WRITE ( 6 , *) '..TIME NORMALIZATION FACTOR. . (TINF) . : ' ,TINF 
READ(4,*) IDATA 

WRITE(6,*) '..0=RELATIVE TEMP; FACTUAL TEMP. . (IDATA) . : ' , IDATA 
READ(4,*) ZI 

WRITE(6,*) '..AXIAL INCREMENT FACTOR. . (ZI) . : ' ,ZI 
READ (4,*) POWER 

WRITE(6,*) '..POWER MULTIPLICATION FACTOR.. (POWER).: ', POWER 
READ(4,*) ACC 
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95 

96 

97 

98 

99 

100 C 

101 C 

102 C 

103 

104 

105 

106 

107 

108 

109 

110 
111 

112 C 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 C 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 
141 

142 C 

143 

144 


WRITE(6,*)' . .OUTPUT ANGLE ACCURACY CRITERIA. . (ACC) ACC 

READ(4,*) TPUMP v , 

WRITE(6,*) '..TYPE OF PUMPING (0=side, l=end)..: , TPUMP ( 

WRITE(6,*) ' , 

WRI TE{ 6 , * ) ' 

...(PROCESSING OF INPUT FILE TWO)... 

WRITE( 6,*) ' ' 

WRITE{ 6 , * ) ' INPUT FILE 12......' 

WRITE(6,*) ' ' 

JJ=0 

5 CONTINUE 

IF((JJ/3)*3.NE.JJ) GOTO 10 
READ{ 5,100, ERR=86) (XX(I) ,1=1,132) 

WRITE( 6 , 100) ( XX ( 1 ) ,1=1,132) 

10 CONTINUE 

IF(JJ.NE.O) GO TO 20 
IF (TPUMP. EQ.O) THEN 
DO 15 J=1 , 80 

READ( 5,100, ERR= 87 ) (XX(1) ,1 = 1, 132) 

WRITE{ G , 100) (XX( I ) ,1=1,132) 

15 CONTINUE 

ELSEIF( TPUMP. EQ. 1) THEN 
DO 16 J= 1 , 81 

READ(5,100,ERR=87) (XX(1) ,1=1,132) 

WRITE (6,100) (XX( I), 1=1, 132) 

16 CONTINUE 
ENDIF 

20 JJ=JJ+1 
P=P+1 

READ(5, 100 ,ERR=88) (XX(I) ,1=1,132) 

WRITE( 6 , *) ' BLOCK OF TEMPERATURE PROFILES FROM INPUT FILE #2 
READ( 5 , 100 ) (XX(I) ,1=1,132) 

WRITE( 6 , 100) (XX ( I ) ,1=1,132) 

READ (5,100) (XX(I) ,1=1,132) 

WRITE (6,100) (XX(I) ,1=1,132) 

READ( 5 ,105) (XX( I ) ,1=1,50) ,DIFF 
WRITE( 6,105) (XX(l) ,1=1,50) ,DIFF 
READ ( 5 , 105) (XX(I) ,1=1,50) , RADI 
WRITE (6,105) (XX(I) ,1=1,50) ,RADI 
READ( 5,105) (XX(I) ,1=1,50), LENG 
WRITE( 6,105) (XX(I) ,1=1,50) ,LENG 
READ( 5 , 105 ) (XX(I) ,1=1,50) , SEC 
WRITE{ 6,105) (XX(l) ,1=1 ,50) , SEC 

DO 25 J=l,2 

READ( 5 , 100 ) ( XX ( I ) ,1=1,132) 
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145 

146 

147 

148 

149 

150 C 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 C 

173 C 

174 C 

175 

176 

177 

178 

179 

180 
181 
182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 


WRITE ( 6 , 100 ) (XX(I) , 1=1, 132) 
25 CONTINUE 


READ(5,115) ( XX ( I ) , 1=1,7) , (LON(I) , 1=1, N) 
WRITE( 6,115) { XX{I ), 1=1,7 ),(LON{I), 1=1, N) 


DO 30 J=1,M 

READ{ 5 , 125 ) (XX(l), 1=1,2) , RAD{ J) , (TEMP{I, J) , 1=1, N) 

WRITE( 6 , 125 ) { XX {I),I=l,2) f RAD { J ) , { TEMP (I,J),I=1,N) 

30 CONTINUE 
DO 40 J=1,M 
DO 35 1=1, N 

TEMP(I, J)=TEMP(I, J)*POWER 
IF{ IDATS . EQ . 0 ) TEMP (I , J ) = TEMP { I , J) *TNF 
35 CONTINUE 
40 CONTINUE 
WRITE{ 6 , * ) ' 

2S3M1! . THERMAL LENSING PROGRAM OUTPUT DATA' 

WRITE( 6, * ) 

WRITE( o'*) ' ' INDEX OF REFRACTION MATRIX ' 

DO 50 J=1,M 
DO 45 1=1, N 

INDEX ( I , J ) = IOR+ ( DT*TEMP ( I , J) ) 

45 CONTINUE 

WR ITE( 6,125) (XX(I) ,1=1,2) ,RAD{J) , (INDEX(I,J) ,I=1,N) 

50 CONTINUE 


...{CALC OF LENSING AND DISTORTION)... 

PI=3 . 141592653589793D00 
AREA=PI*ZI**2 
SIOR0=0 
DO 60 J= 1 , M- 1 

angij=o.odoo 

SIOR J=0 . 0D00 
SINTJ=O.ODOO 
DO 55 1=1, N-l 
IF(I.EQ.IO) GOTO 65 
IF(J.EQ.l) R1=-RAD(J+1) 

IF(J.GT.l) R1=RAD( J-l) 

R2=RAD(J) 

R3=RAD(J+1) 

IF(J.EQ.l) T1=TEMP( I , J+l) 

IF(J.GT.I) T1=TEMP{ I, J-l) 

IF(J.EQ.l) T 4 = TEMP (I+1,J+1) 

IF(J.GT.l) T4=TEMP(I+1, J-l) 

IF(J.EQ.l) T7=TEMP(I+2,JU) 

IF(J.GT.l) T7=TEMP{I+2, J-l) 

T2=TEMP ( I , J ) 
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195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 


T5=TEMP ( I + 1 , J ) 

T8=TEMP(I+2,J) 

T3=TEMP(I , J+l) 

T6=TEMP(I+1, J+l) 

T9=TEMP(I+2,3+l) 

Zl=LON(I) 

Z2=LON(l+l) 

Z3=LON(I+2) 

Z = LON (I)+(ZI/2) 

CALL INTERP(Z,Z1,Z2,Z3,T,T1,T4,T7 ,A,B,C) 

T1C=T 

CALL INTERP{Z,Z1,Z2,Z3,T,T2,T5,T8,A,B,C) 

T2C=T 

CALL INTERP(Z ,Zl,Z2,Z3,T,T3 ,T6 ,T9,A,B,C) 

T3C=T 

Z=L0N(I+l)+(ZI/2) 

CALL INTERP(Z,Zl,Z2,Z3,T,Tl,T4,T7,A,B,C) 

T4C=T 

CALL INTERP(Z,Z1,Z2,Z3,T,T2,T5,T8,A,B,C) 

T5C=T 

CALL INTERP(Z,Z1,Z2,Z3,T,T3,T6,T9,A,B,C) 

T6C=T 

65 CONTINUE 
LONI=LON(I) 

IF( 1 . EQ , 1) THEN 

CALL INITIAL(RO,ANGOJ) 

SUMI OR=N 1 +-N2 * SQRT ( l+( ( ( (R0-RM)*RADI)/(ZI/2*LENG) )**2 ) ) 

SIORJ=SIORJ+SUMIOR 

IF(J.EQ.l) SIORO=SIORJ 

CALL MIDPT(RO,ANGOJ) 

ELSEIF(I.LT.N-l) THEN 

CALL FINAL(RO,ANGOJ) 

SUMI0R=N1*SQRT( l+( ( ( (RM-RI )*RADI)/(ZI/2*LENG) )**2) ) + 

1 N2*SQRT(1+((((R0-RM)*RADI)/(ZI/2*LENG))**2)) 

SIOR3=SIOR3 +SUMIOR 
IF(J.EQ.l) SIORO=SIORJ 
CALL MIDPT(RO, ANGOJ) 

ELSEIF(I.EQ.N-l) THEN 

CALL FINAL(RO, ANGOJ) 

REND{ J)=RO 

SUMIOR=Nl*SQRT(l+(( ( (RM-RI )*RRDI )/( ZI/2*LENG )) **2 ) ) + 
1 N2*SQRT(l+( ( ( (R0-RM)*RADI)/(ZI/2*LENG))**2) J 

SIORJ=SIORJ+SUMIOR 
N1=N2 

IF(J.EQ.l) THEN 
DISTOR(J)=0 
ELSEIF(J.GT.l) THEN 

DISTOR(J)=100*((RO**2-R2**2)/R2**2) 

ANGO( J)=ASIN(N1*SIN(ANG0J) ) 

ENDIF 
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245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

274 

275 

276 

277 

278 

279 

280 
281 
282 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 


IF(J.EQ.l) SIORO=SIORJ 
IF( S IOR J-SIORO . EQ . 0 } THEN 
FOCAL { J ) =0 

ELSE IF{ SIOR J-SIORO . NE . 0 ) THEN 

FOCAL( J)=-(RO*RADI)**2/(ZI*LENG*(SIORJ-SIORO) ) 

END IF 

ENDIF 

55 CONTINUE 
60 CONTINUE 
WRITE(6,*)' 

WRITE( 6 ,*) ' 

WRITE( 6,135) 

WRITE(6,*)' 

WRITE{ 6 , *) ' 

1 

WRITE(6,*) ' 

1 

WRITE(6,*) ' 

1 

DO 70 J=1,M 
WRITE(6,130) 

70 CONTINUE 
WRITE(6,*)' 

AFOCAL=0 
DO 75 J=1,M 
AFOCAL=AFOCAL+FOCAL( J ) 

75 CONTINUE 

IF(AFOCAL.EQ.O) GOTO 76 
AFOCAL= AFOCAL/ (M-2) 

76 CONTINUE 
ADISTOR^O 
DO 80 J=1,M 

IF(DISTOR(J) .GT.ADISTOR) ADISTOR=DISTOR( J) 

80 CONTINUE 

WRITE(6,140) AFOCAL/ 100 
WRITE{ 6,145) ADISTOR 
HRITE(6,*)' ' 

WRITE(6,*)’ END THERMAL LENSING OUTPUT..’ 

WRITE ( 6 , * ) ’ 

HRITE( 6,*) 1 ’ 

I F( AFOCAL. EQ.0) THEN 
IF(P.GT.l) THEN 

WRITE(6, *)’... PROGRAM TERMINATED. LASER ROD HAS ESSENTIALY’, 
1 ’ RELAXED TO STATIC CONDITIONS.’ 

GOTO 90 
ENDIF 

ENDIF 
GOTO 5 
86 CONTINUE 

WRITE(6,*) ’ ERROR ON LINE 109’ 

GOTO 90 


WL,TOI,SEC 

RAY ANGLE OUT RADIUS FOCAL ’ 

’ DISTORTION ’ 

(no.) (radians) (cm) L ENG (cm) ’ 

‘ (%)* 


» t 

RAD( J) , ANGO( J) ,REND( J) ,FOCAL( J) ,DISTOR( J) 
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296 

297 

298 

299 

300 

301 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 

332 

333 

334 

335 

336 

337 

338 

339 

340 

341 

342 

343 

344 

345 


87 CONTINUE 
WRITE (6,*) 
GOTO 90 

88 CONTINUE 
WR1TE{6,*) 


ERROR ON LINE 116/121' 


ERROR ON LINE 127' 


100 F0RMAT(132A1) 

105 FORMAT ( 50 A1 , Ell . 6 ) 

115 FORMAT (7ftl, 14F12 .2) 

125 FORMAT( 2A1 ,F5 . 2 , 14E12 .5) 

130 FORMAT(7X,F5.2,17E15.7) 

135 FORMAT(//,22X, 'THERMAL LENSING CALCUATIONS' , 

1 / , 15X, 'WAVELENGTH OF INTEREST(um) \E12.6, 

1 /,15X, 'TEMPERATURE OF INTEREST(deg C).. \E12.6, 

1 / , 15X, 'TIME (sec) \E12.6) 

140 FORMAT ( 7 X, 'AVERAGE OF THE FOCAL LENGTHS (meters) =',E15.6) 

145 FORMAT (7X, 'INCREMENTAL DISTORTION (%) RANGES FROM 0.0 TO ,E10.3) 


90 CONTINUE 
END 


C 

C 

SUBROUTINE INITIAL(RO,ANGOJ) 


C 

c 

C This subroutine calculates the initial entry into the temperature 
C gradient matix and perdicts the initial angle and radius out for 
C the first axial increment of the laser rod. 

C 

REAL IOR,LONI ,N1,N2 
C 

DOUBLE PRECISION ANGMJ,ANGIJ,ANGOJ,Al,A2,A3,RI ,RM,RO,Z,A,B,C 


COMMON /BLOCKl / IOR , DT , ACC ,N1 , N2 , ZI , RI ,RM ,TM,LONI , I 
COMMON /BLOCK2/R1 ,R2 ,R3,Z1,Z2 ,Z3 
COMMON /BLOCR3/Tl,T2,T3,T4,T5,T6,T7 ,T8,T9 
COMMON /BLOCK4/T1C , T2C , T3C , T4C , T5C , T6C , T10 , Til , T12 

RI =R2 

CALL INTERP(RI,R1,R2,R3,T,T1C,T2C,T3C,A,B,C) 

DTR=(2*A*RI)+B 

Z=LONI +ZI / 2 

CALL INTERP (Z,Zl,Z2,Z3,T, T2 ,T5,T8,A,B,C) 
DTZ=(2*A*Z)+B 
IF(DTR.EQ.O) THEN 
ANGMJ=0,0D00 
ELSEIF(DTR.NE.O) THEN 
ANGMJ=ASIN(1 .0) 

ENDIF 

IF(DTZ . GE . I . 0E- 30 ) ANGMJ=ATAN(DTR/DTZ) 
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346 RM=R2 

347 TM=T2C 

348 Nl=IOR+DT*( (TM+T2 ) / 2 ) 

349 N2=I0R+DT*( (TM+T5)/2) 

350 A1=SI N ( ANGMJ- ANGI J ) 

351 A2=(N'2/N1)*A1 

352 A3=ASIN ( A2 ) 

353 ANGOJ=ANGMJ-A3 

354 C ANGOJ= ANGMJ- ASIN ( (N2/N1)*SIN{ ANGMJ- ANGI J ) ) 

355 RO=(ZI/2)*TAN(ANGOJ)+RM 

356 CALL REFINE(R0,ANGI3, ANGMJ, ANGOJ,T4,T5,T6) 

357 C 

358 RETURN 

359 END 

360 C 

361 C 

362 SUBROUTINE MIDPT(RO,ANGOJ) 

363 C 

364 C This subroutine calculates the intermediate angles and predicts 

365 C the mid-point position of the next axial increment prior to 

366 C entry into the next increment of the temperature gradient matrix. 

367 C 

368 REAL IOR,LONI ,N1 ,N2 

369 C 

370 DOUBLE PRECISION ANGMJ,ANGIJ,ANGOJ,Al,A2,A3,RI,RM,RO,Z,A,B,C 

371 C 

372 COMMON /BLOCKl/IOR,DT,ACC,Nl,N2,ZI,RI,RM,TM,LONI,I 

373 COMMON /BL0CK2/R1,R2,R3,Z1,Z2,Z3 

374 COMMON /BLOCK3/Tl,T2,T3,T4,T5,T6,T7,T8,T9 

375 COMMON /BLOCK4/T1C,T2C,T3C,T4C,T5C,T6C,T10,T11 I T12 

376 C 

377 

378 RM=RO 

379 ANGI J=ANGOJ 

380 CALL INTERP(RM,Rl,R2 ,R3,T,T1 ,T2,T3,A,B,C) 

381 T10=T 

382 CALL INTERP(RM,R1 ,R2,R3,T,T4,T5,T6,A,B,C) 

383 T11=T 

384 DTR=(2*A*RM)+B 

385 CALL INTERP(RM,R1,R2 ,R3,T,T7 ,T8,T9, A,B,C) 

386 T12=T 

387 Z=LONI+ZI 

388 CALL INTERP(Z,Z1,Z2,Z3,T,T10,T11,T12,A,B,C) 

389 DTZ = ( 2*A*Z ) +B 

390 TM=T 

391 IF(DTR.EQ.O) THEN 

392 ANGMJ=0 .0D00 

393 ELSEIF(DTR.NE.O) THEN 

394 ANGMJ=ASIN(1.0) 

395 ENDIF 


74 


396 IF(DTZ.GE.1.0E-30) ANGMJ=ATAN ( DTR/DTZ ) 

397 CALL INTERP(RM,Rl,R2,R3,T,T4C,T5C,T6C,A,B,C) 

398 N1=N2 

399 N2=I0R+DT*{ (TM+T)/2) 

400 A1=SIN( ANGMJ-ANGI 0 ) 

401 A2= ( N2/K1 ) *Al 

402 A3=ASIN(A2) 

403 ANGOJ=ANGMJ-A3 

404 C ANG0J=ANGMJ-ASIN({N2/N1)*SIN(ANGMJ-ANGI J) ) 

405 R0=(ZI/2)*TAN(ANG0J)+RH 

406 CALL REFINE(RO, ANGI J,ANGMJ, ANGOJ ,T4C,T5C,T6C) 

407 C 

408 RETURN 

409 END 

410 C 

411 C 

412 SUBROUTINE FINAL (RO, ANGOJ) 

413 C ■ £ 

414 C This subroutine calculates the final angle and radius out for 

415 C all subsequent axial increments of the laser rod. 

416 C 

417 REAL IOR,LONI,Nl,N2 

418 C 

419 DOUBLE PRECISION ANGMJ,ANGIJ,ANGOJ,Al,A2,A3,RI,RM,RO,Z,A,B,C 

420 C 

421 COMMON /BLOCKl/IOR,DT,ACC,Nl,N2,ZI,RI,RM,TM,LONI,I 

422 COMMON /BL0CK2/R1,R2,R3,Z1,Z2,Z3 

423 COMMON /BLOCK3/Tl,T2,T3,T4,T5,T6,T7,T8,T9 

424 COMMON /BLOCK4/T1C,T2C,T3C,T4C,T5C,T6C,T10 ,T11,T12 

425 C 

426 

427 RI=RM 

428 RM=R0 

429 ANGI J=ANGOJ 

430 IF(I.LT.IO) CALL INTERP(RM,R1 ,R2 ,R3 ,T ,T1C,T2C,T3C,A,B,C) 

431 IF(I.EQ.IO) CALL INTERP(RM R1,R2,R3,T,T4C,T5C,T6C,A,B,C) 

432 DTR=(2*A*RM)+B 

433 CALL INTERP(RM,Rl ,R2 ,R3,T,T1,T2,T3 , A,B,C) 

434 T10 =T 

435 CALL INTER? (RM,Rl I R2,R3,T,T4 I T5,T6 r A,B,C) 

436 T11=T 

437 CALL INTERP(RM,R1,R2 ,R3,T,T7 ,T8,T9,A,B,C) 

438 T12=T 

439 IF(I.LT.IO) Z=LONI+ZI/2 

440 IF(I.EQ.IO) Z=LONI+ZI/2 

441 CALL lNTERP(Z,Zl,Z2,Z3,T,T10 r Tll I T12,A,B f C) 

442 DTZ=(2*A*Z)+B 

443 IF(DTR.EQ.O) THEN 

444 ANGMJ=0.0DQ0 

445 ELSEIF(DTR.NE.O) THEN 
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ANGMJ= ASIN (1.0) 


446 

447 END IF 

448 IF(DTZ .GE. 1.0E-30) ANGMJ= ATAN ( DTR/DTZ ) 

449 TM=T 

450 N1=N2 

451 IF(I.LT.IQ) N2 = IOR+DT* { ( TM+T11 ) / 2 ) 

452 IF(I.EQ.IO) N2= IOR+DT* ( ( TM+T12 )/2 ) 

453 A1=SIN(ANGMJ-ANGIJ) 

454 A2=(N2/N1)*A1 

455 A3=AS IN { A2 ) 

456 ANGOJ=ANGMJ-A3 

457 C ANGOJ=ANGM J-ASIN( (N2/N1 ) *SIN( ANGMJ-ANGI J) ) 

458 RO=(ZI/2)*TAN(ANGOJ)+RM 

459 IF(I.LT.IO) CALL REFINE(R0,ANGIJ,ANGMJ,ANG0J,T4,T5,T6) 

460 IF(I.BQ.IO) CALL REFINE ( RO , ANGI J , ANGM J , ANGO J , T7 , T8 , T9 ) 

461 C 

462 RETURN 

463 END 

464 C 

465 C 

466 SUBROUTINE INTERP(R,R1,R2,R3,T,T1,T2,T3,A,B,C) 

467 C 

468 C 

469 C The subroutine INTERP is a binomial interpolation routine using 

470 C any three points in the R (radial) or Z (axial) direction to find 

471 C the coefficents of A,B, and C of the equation AX**2 + BX 4 C = T 

472 C so that a temperature T can be determined for any X (R or Z) 

473 C within the non-linear temperature matrix. 

474 C 

475 DOUBLE PRECISION R,A,B,C 

476 C 

477 A= ( ( T1 -T2 ) * (R1-R3 ) - ( T1 -T3 ) *(R1-R2 ) ) / ( (R1-R3 )*(R1**2-R2**2 )— 

478 1 (R1-R2)*<R1**2-R3**2) ) 

479 B={ (Tl-T3)*(Rl**2-R2**2)-{Tl-T2)*(Rl**2-R3**2))/( (R1-R3)* 

480 1 (R1**2-R2**2)-(R1-R2)*(R1**2-R3**2) ) 

481 C=T1-A*R1**2-B*R1 

482 T=A*R**2+B*R+C 

483 C 

484 RETURN 

485 END 

486 C 

487 C 

488 SUBROUTINE REFINE(R0, ANGI J, ANGM J,ANG0J,T1,T2,T3) 

489 C 

490 C 

491 C The subroutine REFINE is a continous numerical iteration routine 

492 C used to recalculate the final thermal lensing output angle for 

493 C each specified axial increment within each ray until the user 

494 C specified accuracy criteria for the output angle are met. 

495 C 


76 



496 C 

497 REAL I0R,L0NI,N1,N2 

498 C 

499 DOUBLE PRECISION ANGMJ,ANGI J,ANGOJ,Al,A2,A3,RI,RM,RO 

500 DOUBLE PRECISION ROO,D,A,B,C 

501 C 

502 COMMON /BLOCK1/IOR,DT,ACC,N1,N2,ZI,RI,RM,TM,LONI,I 

503 COMMON /BL0CK2/R1,R2,R3,Z1,Z2,Z3 

504 C 

505 L=0 

506 5 CONTINUE 

507 L=L + 1 

508 CALL INTERP(R0,R1,R2,R3,T,T1,T2,T3,A,B,C) 

509 N2=I0R+DT*((TM+T)/2) 

510 A1=SIN(ANGMJ-ANGIJ) 

511 A2-(N2/N1)*A1 

512 A3=ASIN(A2) 

513 ANGOJ=ANGMJ-A3 

514 C ANG0J=ANGMJ-ASIN((N2/N1)*SIN(ANGMJ-ANGIJ)) 

515 R00= { Z I / 2 ) *TAN { ANGO J ) +RM 

516 D=RO 

517 R0=R00 

518 IF( ABS(D-ROO) .GT.ACC) GOTO 5 

519 C 

520 RETURN 

521 END 
No errors detected 
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Appendix B 


Program Input Files 


1. Thermal Lensing Model Input File #1 

This file (THERLEN 1 . I NP ) contains thermal lensing 
material char acter i st i cs and other constants required to 
control the Thermal Lensing Model. The example data 
depicted below are for a pulsed side-pumped Nd : YAG rod. 


File 


THERMAL LENSING TNPUT DATA FOR Nd : YAG ROD 


1.15 J WAVELENGTH OF INTEREST (MICRONS) 

37.5 ; TEMPERATURE OF INTEREST (24-42 DEGREES C) 

1.825 ; INDEX OF REFRACTION 

8.90E-06 ; ALPHA (DN/DT) 

0.33 ; TEMPERATURE NORMALIZATION FACTOR 

0.05 ; TIME NORMALIZATION FACTOR 

1 ; 0 - RELATIVE VALUES IN INPUT FILE; 1 - ACTUAL VALUES 

0.10 : AXIAL INCREMENT 

1.0 ; POWER MULTICATTON FACTOR 

1.0E-06 ; OUTPUT ANGLE ACCURACY CRITERIA 

0 ; TYPE OF PUMPING (0=side, l=end) 


2. Thermal Lensing Model Input File #2 

This file ( THERLEN2 . I NP ) is tha. output file from the 
model titled: Temperature Distribution Model for finite 

optically pumped laser rods. It contains that model's 
control criteria and output temperature data blocks versus 
t ime . 

Those data blocks of temperatures contained within 
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this input file are used in the Thermal Lensing Model to 
calculate the thermal lensing and pump beam distortion 
effects. The example data depicted below are for a pulsed 
side-pumped N d : Y AG laser at its steady-state equilibrium 
temperature. Blocks of temperature data are shown from 
3 . 4-3 . 6 seconds . 

I nput File #2 


Input NdtYAG Side _ Pump . . MP(20) . . GA-20 . .GL-None . .Run 001 03/01/90 

MIN. MAX. AND INCREMENT OF RADIUS OF ROD 2.500000E-01 

2.500000E-01 2.500000E-01 ... , n / nnnnn 

MIN , MAX AND INCR. OF ASPECT RATIO (L/A) 30.400000 

30.400000 50.000000 

INVERSE OF RADIAL THERMAL DIFFUSIVITY 26.947000 

THERMAL CONDUCTIVITY 1.030000E-01 1.030000E-01 

MIN. RADIAL COOLING PARAMETER 20.000000 

MAX. RADIAL COOLING PARAMETER 20.000000 

TYPE OF SCALE 1(L0G) OR 0( LINEAR) nnnnnn 

*™ . l NC 5f RATIO^OR^VALUES 6 . 500000E-04 

6 . 500000E-04 , 

TYPE OF SCALE l(L0G; OR 0( LINEAR) 1 

. . Atro .l NC §fS EN U S K°M0i F CmLg RATIO 6.500000E-04 

6 . 56 OOOOE-O 4 , 

TYPE OF SCALE l(LOG) OR 0( LINEAR) 1 

AND INCREMENT (SET TO 1 IF LOG) 1.000000 

. . . TYPE OF DATA RATI0=0 ACTUAL=.l 
PARAM=2 

PROPORTIONAL TO L/A=3 2 

VARIATION OF END DATA SEPERATE=0 

T0GETHER=1 1 

!!!!!' MAxiMUM NUMBER OF ROOTS . . . RADIAL 17 

MAXIMUM NUMBER OF ROOTS . . . L0NGIT 17 

NUMBER OF ORTHOGONAL CHECKS RADIAL 2 

NUMBER OF ORTHOGONAL CHECKS L0NGIT 3 

BESSEL ROOTS SCHEME ... 1 COMPUTED 

2 GIVEN 

3 ALL (GIVEN <21) 

(COMPUTED GE 21) 
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datntat R ?P TS SCHEME is 9 

^TYpE 0? N ^ N SS ATION “VISIONS 

radial eigenfunction 

SQUARE PULSE 
GAUSSIAN 
rA , TrTT PARABOLIC 

LONGIT EIGNFUNCTION 

CONSTANT 

CONSTANT & ATTENUATED +V 

focused(defocused) +2 

^RPPSSP^^.^TENUATED +3 


1000 


-ORDER 

+0 

+1 

+2 

-ORDER 

+0 


NOTE 


^Pr^RS^OSEPFECKALLOWED 


PUMP BEAM €0DE SIAN ENVELO g E ASSUMED 

:.opi™res N DATA - 


0 

1 

2 


NUMBER OF RADIAL PROFILES 

NUMBER OF LONGIT PROFILES 

P?E^ 0US PUMPPING 

SINGLE PULSE 

OR MULTI-PULSED 

;;fU • REL. TIME OF INTEREST' 

iPRE-O^A?? L iJ^ 0F ^EREST 

AmS!^® ?§ ? C ^ N gggJ 

MULTIPLE PULSES G; 2 

ONTIME PER PULSE 1 ioOOOOF-OA 

OFFTIME PER p ULSE 4 ' 987000F-0? 

ONTIME SAMPLING RATE H/000E 02 

;a ° p ™ E sampling rate 

ABSORBED ENERGY PER PULSE J/CM3 

DEBUG DESIRED 7 

mtm ’ ' m ’ av ' NOT DESIRED 

• .‘min; max ’ longitudinal^^ • 000000E+00 

■•RADIAL REFERENCE POINTS LCULATI0N RANGE . . 0.000000E+00 

■ .AXIAL REFERENCE POINTS * " ' -0 0.000000E+00 

0 O.OOOOOOE+OO 

ROD ASPECT RATIO (L/A) 

- ■ - 30.400000 

************************** *********r*rr 

PARAMETER GL OF FRONT 6.500000E-04 

M5!L?L?LlftT,„„| ii |oooooE-o4 ie 


11 

11 


3.400000 
3 . 600000 
1 

1.000000 
2 


1 
1 

1574.280000 
0, NDEBUG= 


100 


1 

1 . 000000 
1.000000 
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TEMPERATURE DISTRIBUTION 


RAO I Al 0 1 FfUSl V1TY 2694706+02 

RADIUS Of ROO 250000 

LENGTH OF ROO 7.600000 

TIME IN SECONDS 345013E+01 


L/R 

.00 

.10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

.90 

1.00 

.00 

.723826+00 

.72383E+00 

.72383E+00 

.72383E+00 

.723836+00 

, 72383E+00 

.72383E+00 

.723836+00 

.723836+00 

.72383E+00 

.723826+00 

.10 

.713306+00 

.71331E+00 

.71331E+00 

.713316+00 

.713316+00 

.71331E+00 

.71331E+00 

.713316+00 

.71331E+00 

.71331E+00 

.71330E+00 

.20 

.694036+00 

.694046+00 

.694046+00 

.69404E+00 

.694046+00 

.69404E+00 

.694046+00 

.69404E+00 

.69404E+00 

.69404E+00 

.69403E+00 

.30 

.663206+00 

.66320E+00 

.663206+00 

.663206+00 

.6632QE+00 

.66320E+00 

.663206+00 

.66320E+00 

.663206+00 

.66320E+00 

.663206+00 

.40 

.62046E+00 

.62046E+00 

.620466+00 

.620466+00 

.620466+00 

.62047E+00 

.620466+00 

.62046E+00 

.620466+00 

.620466+00 

.620466+00 

.50 

.565216+00 

.56522E+00 

.565226+00 

.56522E+00 

.565226+00 

.56522E+00 

.565226+00 

.565226+00 

. 56522E+00 

.56522E+00 

.56521E+00 

.60 

.496886+00 

. 49689E+00 

.496896+00 

.49689E+00 

.49689E+00 

.496896+00 

.496896+00 

.496896+00 

.496896+00 

.496896+00 

.496886+00 

.70 

.415526+00 

.41552E+00 

.415526+00 

.41552E+00 

.41552E+00 

.415526+00 

.415526+00 

.41552E+00 

.415526+00 

.41552E+00 

.415526+00 

.80 

.322886+00 

.32288E+00 

•32288E+00 

.322886+00 

.322586+00 

.322886+00 

.322886+00 

.322886+00 

.322886+00 

.322886+00 

.322886+00 

.90 

,221906+00 

.221906+00 

.221906+00 

.221906+00 

.221906+00 

.221906+00 

.22190E+00 

.22190E+00 

.221906+00 

.22190E+00 

.221906+00 

.00 

.10556E+00 

. 105566+00 

.105566+00 

. 10556E+00 

. 1 0556E+00 

. 105566+00 

.105566+00 

.105566+00 

.105566+00 

. 1 0556E+00 

. 10556E+00 


TEMPERATURE DISTRIBUTION 


RADIAL D I F FUSIVITY 2694706+02 

RADIUS OF ROO 250000 

LENGTH OF ROO 7.600000 

TIME IN SECONDS 350000E+01 


L/R 

.00 

.10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

.90 

1.00 

.00 

.646276+00 

64627E+00 

.646276+00 

.646276+00 

.64627E+00 

.64627E+00 

.64627E+00 

.646276+00 

.646276+00 

.64627E+00 

.64627E+00 

.TO 

.640066+00 

.640Q7E+00 

.64007E+00 

.64007E+00 

.64007E+00 

.64007E+00 

.64007E+00 

.64007E+00 

.64007E+00 

.64007E+00 

.640066+00 

.20 

.621446+00 

.62145E+00 

.621456+00 

.621456+00 

.621456+00 

. 62145E+00 

.621456+00 

.621456+00 

.62145E+00 

.62*456+00 

.621446+00 

.30 

.590416+00 

.590416+00 

.590416+00 

. 59041E+00 

.590416+00 

.59041E+00 

.59041E+00 

.590416+00 

.590416+00 

.590416+00 

.590416+00 

.40 

.546936+00 

.54694E+00 

.546946+00 

.546946+00 

.54694E+00 

.54694G+00 

.54694E+00 

.546946+00 

.54694E+00 

.546946+00 

.54693E+Q0 

.50 

.49096E+00 

.490966+00 

.490966+00 

.490966+00 

.490966+00 

.49096E+00 

.490966+00 

.490966+00 

.490966+00 

.490966+00 

.490966+00 

.60 

.422446+00 

.42244E+00 

.42244E+0Q 

.42244E+00 

.422446+00 

.42244E+00 

.422446+00 

.422446+00 

.422446+00 

.422446+00 

.42244E+00 

.70 

.341696+00 

.34169E+00 

.341696+00 

.341696*00 

.34169E+0Q 

.34169E+00 

.341696+00 

.341696+00 

.341696+00 

.341696+00 

.341696+00 

.80 

.250086+00 

.250086+00 

.250086+00 

.250086+00 

.250086+00 

.25008E+00 

.250086+00 

.25008E+00 

.25008E+DG 

.250086+00 

.250086+00 

.90 

.150936+00 

.150936+00 

.15093E+00 

.150936+00 

.150936+00 

. 15093E+00 

.150936+00 

.150936+00 

.150936+00 

.150936+00 

.150936+00 

1.00 

.498566-01 

.498566-01 

.498566*01 

.49857E-01 

.498576*01 

.49857E-01 

.49857E-Q1 

.498576*01 

.498566-01 

.498566-01 

.498566-01 


TEMPERATURE DISTRIBUTION 


RADIAL D I FRJSI VITY 2694706+02 

RADIUS OF ROO 250000 

LENGTH OF ROO 7.600000 

TIME IN SECONDS 350013E+01 


L/R 

.00 

.10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

.90 

1.00 

.00 

.723826+00 

.723836+00 

.723836+00 

.723836+00 

.72383E+00 

.723836+00 

.72383E+00 

.723836+00 

.723836+00 

.723836+00 

.723826+00 

.10 

.713306+00 

. 71331 E+OQ 

.713316+00 

.713316+00 

. 71331 E+00 

.71331E+00 

.713316+00 

.713316+00 

.71331 E+00 

.713316+00 

.713306+00 

.20 

.694036+00 

.694046+00 

.694046+00 

.694046+00 

.694046+00 

.694046+00 

.694G4E+00 

.694046+00 

.69404E+0G 

.694046+00 

.694036+00 

.30 

.663206+00 

.663216+00 

.663216+00 

.663216+00 

.663216+00 

.663216+00 

.66321 E+00 

.66321 E+00 

.663216+00 

.663216+00 

.663206+00 

.40 

.620466+00 

.62047E+00 

.620476+00 

.620476+00 

.620476+00 

.620476+00 

.620476+00 

.620476+00 

.620476+00 

.620476+00 

.620466+00 

.50 

.565216+00 

.565226+00 

.56522E+00 

.565226+00 

.565226+00 

.565226+00 

.565226+00 

.565226+00 

.565226+00 

.565226+00 

.565216+00 

.60 

.496886+00 

, 49689E+00 

.496896+00 

.49689E+00 

.496896+00 

.496896+00 

.49689E+00 

.496896+00 

.4*6896+00 

.496896+00 

.496886+00 

.70 

.415526+00 

.41552E+00 

.41552E+00 

.41552E+00 

.415526+00 

.415526+00 

.41552E+00 

.415526+00 

.415526+00 

.415526+00 

.415526+00 

.80 

.322886+00 

.32258E+00 

.322886+00 

.322886+00 

.322886+00 

.322886+00 

.32288E+00 

.322886+00 

.3228SE+00 

.322886+00 

.322886+00 

.90 

.221906+00 

. 22 190E+00 

.221906+00 

.221906+00 

.321906+00 

.221906+00 

.22190E+00 

.221906+00 

.221906+00 

.221906+00 

.221906+00 

.00 

. 105566+00 

. 10556E+00 

.105566+00 

.105566+00 

.105566+00 

.105566+00 

. 1C5566+00 

. 105566+00 

.105566+00 

.105566+00 

.105566+00 
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TEMPERATURE DISTRIBUTION 

RADIAL D I FFUSIVITY 269470E+02 

RADIUS OF ROO 250000 

LENGTH OF RCO 7.600000 

TIME IN SECONDS 3550006+01 


L/R 

.00 

.10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

.90 

1.00 

.00 

.64627E+DQ 

.646276+00 

.64627E+00 

.646276+00 

.646286+00 

.64628E+00 

.64628E+00 

.64627E+00 

.64627E+00 

.646276+00 

.&4627E+00 

.10 

•64Q06E+00 

.64007E+00 

.64007E+00 

.640076+00 

.64007E+00 

.64007E+00 

.64007E+00 

.64007E+00 

.64007E+00 

.64007E+00 

.640066+00 

.20 

.62144E+00 

.62145E+00 

.62145E+00 

.62145E+00 

.62145E+00 

.621456+00 

.62145E+00 

.62145E+00 

.62145E+00 

•62145E+00 

.621446+00 

.30 

.59041E+00 

.59041E+00 

.59041E+00 

.59041E+00 

.59042E+00 

.590426+00 

.59042E+00 

.5904 IE +00 

.590416+00 

.59041E+00 

.59041E+00 

.40 

.54694E+QQ 

.54694E+00 

.54694E+00 

. 54694E+00 

.54694E+G0 

.54694E+00 

.54694E+00 

.5 4694 E+00 

.546946+00 

.54694E+00 

.546946+00 

.50 

.49O96E+0Q 

,49097E+00 

.49097E+00 

.49097E+00 

.49O97E+D0 

.490976+00 

.49097E+00 

.49097E+00 

.490976+00 

.490976+00 

.490966+00 

.60 

.42244E+00 

.42244E+00 

.42244E+00 

.42244E+00 

.42244E+00 

.42244E+00 

.422446+00 

.42244E+00 

.422446+00 

.422446+00 

.422446+00 

.70 

.341696+00 

. 34 169E+00 

.34169E+00 

.34169E+00 

.341696+00 

.341696+00 

.341696+00 

.341696+00 

.34169E+00 

.341696+00 

.34169E+00 

.80 

.25008E+00 

.25006E+00 

.25008E+00 

.25008E+00 

,250086+00 

. 25008E+00 

.250066+00 

.25008E+00 

.25005E+00 

.250086+00 

.250086+00 

.90 

. 15093E+Q0 

. 15093E+00 

.150936+00 

. 15093E+00 

.15093E+00 

. 15093E+00 

.15093E+00 

. 15093E+00 

.150936+00 

. 15093E+00 

. 15093E+D0 

1.00 

.49856E-01 

.49S57E-01 

.49857E-01 

.49857E-01 

.49S57E-Q1 

.498576-01 

.49857E-01 

.49857E-0I 

.49857E-Q1 

.49857E-01 

.498566-01 


TEMPERATURE DISTRIBUTION 

RADIAL DI FFUSIVITY 2694706+02 

RADIUS OF ROO 250000 

LENGTH OF ROO 7.600000 

T IM6 IN SECONDS .3550136+01 


L/R 

.00 

.10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

.90 

1.00 

.00 

.72382E+00 

. 72383E+00 

.72333E+00 

.723836+00 

.72383E+00 

.723836+00 

.723836+00 

.723836+00 

.72383E+00 

.72383E+00 

.72382E+00 

.10 

.713316+00 

.713316+00 

.71331E+00 

.713316+00 

,713516+00 

.71 331 £+00 

.713316+00 

.713316+00 

- 7133 1 E+00 

.71331E+00 

.713316+00 

.20 

.694036+00 

.69404E+00 

.694046+00 

.69404£h 00 

. 69404E+00 

.694046+00 

. 69404 E+00 

.69404E+00 

.69404E+00 

.694046+00 

.69403E+0O 

.30 

.66320E+00 

.663216+00 

•66321E+00 

.663216+00 

.66321E+00 

.66321E+00 

.663216+00 

.66321 E+00 

.6632 1E+00 

.663216+00 

.663206+00 

.40 

. 620466+00 

.62047E+00 

.620476+00 

.620476+00 

.62047E+00 

. 6204 7E+00 

.62047E+00 

.620476+00 

.620476+00 

.620476+00 

.620466+00 

.50 

. 5652 1 E+00 

.565226+00 

. 565226+00 

.565226+00 

.565226+00 

, 56522E+00 

.56522E+00 

-56522E+00 

.56522E+00 

.565226+00 

.56521E+00 

.60 

. 49658E+Q0 

. 49689E+00 

.496896+00 

. 496896+00 

.49659E+00 

. 49689E+00 

.49689E+00 

.496896+00 

.49689E+00 

.496896+00 

.49689E+00 

.70 

.41552E+00 

.41552E+00 

.A1552E+O0 

.41552E+00 

.41552E+00 

.41552E+00 

.415526+00 

.415526+00 

.41552E+00 

.415526+00 

.41552E+00 

.50 

.322886+00 

.32258E+00 

.322886+00 

.322586+00 

.32288E+00 

.32258E+00 

.32288E+00 

.322886+00 

.32288E+00 

.322886+00 

.322886+00 

.90 

.2219GE+00 

.221906+00 

.22190E+00 

.221906+00 

.22190E+00 

. 22190E+00 

.221906+00 

.221906+00 

. 2Z190E+D0 

,221906+00 

.22190E+00 

.00 

.105566+00 

. 10556E+00 

. 10556E+D0 

. 10556E+00 

.105566+00 

. 105566+00 

.105566+00 

. 10556E+D0 

.105566+00 

.10556E+CC 

. 10556E+00 


TEMPERATURE DISTRIBUTION 


RADIAL DI FFUSIVITY 269470E+02 

RADIUS OF ROO 250000 

LENGTH OF ROO 7.600000 

TIME IN SECONDS 3600006+01 


L/R 

.00 

.10 

.20 

.30 

.40 

.00 

.646136+00 

.64613E+00 

.64613E+30 

.646136+00 

.64613E+D0 

.10 

.63993E+0Q 

. 63993E+00 

.63994E+00 

.63994E+00 

.63994 E+00 

.20 

.62131E+Q0 

.621326+00 

.621326+00 

.621326+00 

.621326+00 

.30 

.59C28E+00 

.590286+00 

.59028E+C0 

.590286+00 

.590286+00 

.40 

.546806+00 

.546816+00 

. 546S1E+00 

.546816+00 

.546816+00 

.50 

.490826+00 

.490836+00 

.490836+00 

.49083E+00 

.490836+00 

.60 

.422306+00 

.422316+00 

.42231E+00 

.4223 1 E+00 

.422316+00 

.70 

,341556+OC 

.34156E-HI0 

.341566+00 

.341566+00 

.34156E+00 

.80 

.24995E+00 

.249956+00 

.249956+00 

.24995E+00 

.24995E+00 

.90 

. 15080E+00 

. 15080E+00 

.150806+00 

. 15080E+00 

.150806+00 

.00 

.49752E-01 

.497536-01 

.497536-01 

.497536-01 

.49753E -01 


.50 

.60 

.70 

.80 

.90 

1.00 

.646136+00 

.646136+00 

.64613E+00 

.64613E+00 

.646136+00 

.646136+00 

.63994E+00 

.639946+00 

.63994E+00 

.63994E+00 

.63994E+00 

.639936+00 

.621326+00 

.621326+00 

.62132E+00 

.621326+00 

.62132E+00 

.621316+00 

.59028E+00 

•59028E+00 

.590286+00 

.590286+00 

.59028E+00 

.59028E+00 

.54681E+00 

.54681E+00 

.546816+00 

.546816+00 

.54651E+00 

.546806+00 

.49083E+00 

•49083E+00 

.490836+00 

.490836+00 

. 49083E+00 

.49082E+00 

.42231E-0Q 

.42231 E+00 

.422316+00 

.42231E+00 

.422316+00 

.42230E+00 

.341566+00 

.34156E+00 

.341566+00 

.34156E+00 

.341566+00 

.341556+00 

.249956+00 

.24995E+00 

.24995E+00 

.24995 E+00 

.24995E+00 

.249956+00 

. 15080E+00 

. 15080E+00 

.150806+00 

.15080E+00 

.150806+00 

.150806+00 

.49753E-01 

.49753E-01 

.497536-01 

.49753E-01 

.497536-01 

.497526-01 
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Appendix C 


Program Output File 

Thermal Lensing Model Output File 

The Thermal Lensing Model output file (THERLENS . OUT) 
rewrites the data from each of the input files. It then 
calculates thermal lensing and pump beam distortion effects 
for each block of temperature data obtained from the Tem- 
perature Distribution Model. Currently, the lensing model 
is programmed to calculate 10 equally spaced rays tran- 
siting through the laser material. The output data from the 
Thermal Lensing Model depicted below are for a pulsed side- 
pumped Nd : YAG laser using the input files in Appendix B. 

Thermal Lensing Model Output File 


INPUT FILE in 

THERMAL LENSING INPUT' DATA FOR Nd : YAG ROD 


..WAVELENGTH OF INTEREST IN UM..(WL).I „ 1.150000 

..TEMPERATURE OF INTEREST IN DEG C..(T0l).: 37.500000 

..INDEX OF REFRACTION. . (IOR) . : 1.825000 

. .ALPHA (dn/dt) . . (DT) . : 8 . 900000R-06 

..TEMP NORMALIZATION FACTOR .. (TNF) . : 3.300000E-01 

. .TIME NORMALIZATION FACTOR . . (TINF) . : 5 .OOOOOOE-02 

. . O^RELATIVF, TEMP; 1= ACTUAL TEMP . . (IDATA) . : 1 

..AXIAL I NCREMENT FACTOR . . ( Z I ) . : 1 . OOOOOOE-O 1 

. . POWER MULTIPLICATION FACTOR . . (F0WER) . : 1 . 000000 

..OUTPUT ANGLE ACCURACY CRITERIA . (ACC) . : 1.0000G0E-06 

..TYPE OF PUMPING (0=side, l. = end) . . : O.OOOOOOE+OO 


original 

0f poor 


PAGE is 
QUALITY 
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INPUT FILE # 2 


Input Nd : YAG Side-Pump. .MP(20) . . GA=20. .GL=None. .Run=001 03/01/90 

MIN., MAX. AND INCREMENT OF RADIUS OF ROD 2.500000E-01 

2.500000E-01 2.500000E-01 

MIN . MAX AND INCR. OF ASPECT RATIO (L/A) 30.400000 

30 .400000 50.000000 

INVERSE OF RADIAL THERMAL DIFFUSIVITY 26.947000 

THERMAL CONDUCTIVITY 1.030000E-01 1.030000E-01 

MIN. RADIAL COOLING PARAMETER 20.000000 

MAX. RADIAL COOLING PARAMETER 20.000000 

TYPE OF SCALE l(LOG) OR O(LINEAR) 1 

AND INCREMENT (SET TO 1 . IF LOG) 1.000000 

MIN. AND MAX. FRONT COOLING RATIO OR VALUES 6 . 500000E-04 

6.500000E-04 

TYPE OF SCALE l(LOG) OR 0( LINEAR) 1 

AND INCREMENT (SET TO 1 IF LOG) 1.000000 

MIN. AND MAX. BACK COOLING RATIO OR VALUES 6.500000E-04 

6.500000E-04 

TYPE OF SCALE l(LOG) OR 0( LINEAR) 1 

AND INCREMENT (SET TO 1 IF LOG) 1.000000 

TYPE OF DATA RATIO=0 ACTUAL=1 

PARAM=2 

PROPORTIONAL TO L/A=3 2 

VARIATION OF END DATA SEPERATE=0 

TOGETHER- 1 1 

..... MAXIMUM NUMBER OF ROOTS . . . RADIAL 17 

MAXIMUM NUMBER OF ROOTS . . . LONGIT 17 

NUMBER OF ORTHOGONAL CHECKS RADIAL 2 

NUMBER OF ORTHOGONAL CHECKS LONGIT 3 

BESSEL ROOTS SCHEME ... 1 COMPUTED 

2 GIVEN 

3 ALL (GIVEN <21 ) 

(COMPUTED GE 21) 

ROOTS SCHEME IS - 2 

RADIAL AND LONG. INTEGRATION DIVISIONS 1000 100 

TYPE OF PUMP BEAM 

RADIAL EIGENFUNCTION -ORDER 

SQUARE PULSE +0 

GAUSSIAN + 1 

PARABOLIC +2 
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LONGIT EIGNFUNCTION -ORDER 

CONSTANT +0 

CONSTANT & ATTENUATED +.1 
FOCUSED (DEFCCUSED) +2 

FOCUSED & ATTENUATED +3 
NOTE NO DIFFRACTION EFFECTS ALLOWED 
GAUSSIAN ENVELOPE ASSUMED 
PUMP BEAM CODE 0 0 

ATTENUATION DATA. .0=ATTEN. . 1=ATTEN*L 1 

OPTICAL DATA 3.000000 

NUMBER OF RADIAL PROFILES ... 11 


NUMBER OF LONGIT PROFILES . . . 

CONTINUOUS PUMPPING ... 0 

SINGLE PULSE ... 1 

OR MULTI -PULSED ... 2 

MIN. REL. TIME OF INTEREST 

MAX. REL. TIME OF INTEREST 

TYPE OF SCALE l(LOG) OR 0( LINEAR) 

AND INCREMENT (SET TO 1 IF LOG) 

MULTIPLE PULSES ... 2 

ONTIME PER PULSE 1.300000E-04 

OFFTIME PER PULSE 4.987000E-02 

ONTIME SAMPLING RATE 

OFFTIME SAMPLING RATE 

ABSORBED ENERGY PER PULSE J/CM3 

DEBUG DESIRED ... 1 

NOT DESIRED ... 0. NDEBUG= 

•MIN, MAX, RADIAL CALCULATION RANGE . . O.OOOOOOE+OO 
.MIN, MAX LONGITUDINAL CALCULATION RANGE.. 0.000000E+00 
.RADIAL REFERENCE POINTS... 0 0.000000E+00 

• AXIAL REFERENCE POINTS... 0 0.000000F.+00 


11 


3.400000 
3 . 600000 
1 

1.000000 


1 

1 

1574.280000 


ROD ASPECT RATIO (L/A) 


30.400000 


* * ************ * * * * * * * ******* * * * * * * * * * * * 
PARAMETER GA IS 20.000000 

*************************************** 

PARAMETER GL OF FRONT 6 . 500000E-C4 

PARAMETER GL OF BACK 6 . 500000E-04 

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 


1.000000 

1.000000 
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BLOCK OF TEMPERATURE PROFILES FROM INPUT FILE #2 


TEMPERATURE DISTRIBUTION 

RADIAL DIFFUSIVITY 26947GE+02 

RADIUS OF RCO . 250000E+00 

LENGTH OF ROO 760Q00E+C1 

TIME IN SECONDS 345000E+01 


L/R 

.00 

.10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

.90 

1.00 

.00 

.64626E+00 

.646276+00 

.646276+00 

.64627E+00 

.64627E+00 

.64627E+00 

.64627E+00 

.64627E+00 

.64627E+00 

.64627E+00 

.646266+00 

.10 

.640066+00 

.640Q6E+QG 

.64006E+00 

.640056+00 

.64006E+00 

.64006E+00 

.640066+00 

.640066+00 

.640066+00 

.64006E+00 

.640066+00 

.20 

.62K4E+00 

.62144E+Q0 

.621446+00 

.62144E+00 

.62144E-00 

.621446+00 

.62144E+00 

.621446+00 

.621446-00 

.62144E+00 

.621446-00 

.30 

.590406+00 

.59041 £+00 

.59041E+00 

.59041E+00 

.59041E+00 

.5904 IE-00 

.590/1E+00 

.59041E+00 

.590416+00 

.59041E+00 

.590406+00 

.40 

.546936+00 

.546946+00 

.546946+00 

.54694E+00 

.54694E+00 

,546946+00 

.546946+00 

.54694E+00 

.54694E+Q0 

.54694E+00 

.546936+00 

.50 

.490966+30 

.49096E+00 

.49096E+00 

.49096E+00 

.49096E+00 

.49096E+00 

.49096E+00 

.49096E+00 

.490966*00 

,49C96E*-00 

.490966+00 

.60 

.42244E+D0 

.42244E+00 

.422446+00 

.42244E+00 

.42244E+00 

.422446+00 

.42244E+00 

.422446+00 

.422446+00 

.422446+00 

.42244E+00 

.70 

.341696+00 

.341696+00 

.341696+00 

.34169E+00 

.341696+00 

.341696+00 

.341696+00 

.34169E+00 

.341696+00 

,341696+00 

.341696+00 

.80 

. 250086+00 

.25008E+00 

. 25008E+00 

.250066+00 

,250036-00 

.250086+00 

.25Q08E+00 

.25008E+00 

.250086+00 

.250086+00 

.250086-00 

.90 

. 150936+00 

.15093E+00 

.150936+00 

.150936+00 

.15093E+00 

. 15093E+00 

.150936+00 

. 15093E+00 

. 150936 *-00 

.150936+00 

.150936-00 

.00 

. 49856E-01 

.498566-01 

.498566-01 

.49856E-01 

.498566-01 

.498566-01 

.498566-01 

.49856E-01 

.49856E-01 

.498566-01 

.498566-01 


THERMAL LEHSING PROGRAM OUTPUT DATA 


INDEX OF REFRACTION MATRIX 


.00 

.18250E+01 

.182506+01 

.182506+01 

.18250E+01 

. 13250E+01 

. 18250E+01 

.182506+01 

.182506+01 

.182506+01 

. 18250E-C1 

.182506+01 

.10 

.182506-01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

.18250E+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

20 

. I8250E+01 

.182506+01 

.182506+01 

.182506+01 

.18250E+01 

.182506-01 

.132506+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

.30 

. 15250E+G1 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

. 18250E+01 

.182506+01 

. 18250E+01 

.182506+01 

.18250E+01 

.18250E+01 

.40 

.182506+01 

. 13250E+01 

.182506+01 

.182506+01 

. 18250E+01 

. 18250E+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

. 18250E+Q1 

.50 

.182506-01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

.60 

.182506+01 

.132506+01 

.182506+01 

. 182596+01 

.182506+01 

. 18250E+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

.70 

. 162506+01 

.18250E+01 

.182506+01 

. 18250E+01 

.182506+01 

.18250E+01 

.182506+01 

,182506+01 

.182506+01 

,182506+01 

.132506+01 

.80 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

182506+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

. 18250E+01 

.90 

. 162506-01 

.182506+01 

. 18250E+01 

.',a25CE+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

. 18250E+01 

.182506+01 

.182506+01 

,00 

. 18250E+01 

. 182506+01 

.182506+01 

. 18250E+01 

.182506+01 

.182506+01 

.182506+01 

. 1825CE+01 

.182506+01 

.182506+01 

. 18250E+01 


THERMAL LENS I NG CALCUATIONS 

WAVELENGTH OF INTEREST(um) 1150G06+01 

TEMPERATURE OF INTEREST (deg C).. .375000E+02 

TIKE (see) 3450006+01 


RAY 

ANGLE OUT 

RADIUS 

fcx:al 

DISTORTION 

(no. ) 

(radians) 

(cm) 

L ENG 'em) 

(X) 

.00 

. OOOOOOQE* , CQ 

.OOOOOOOE+OO 

.OOOOOOOE+OO 

.00000006+00 

.10 

.00000006+00 

.10000006+00 

.OOOOOOOE+OO 

, OOOOOOOE+OO 

.20 

.00000006+00 

. 200QQ00E+00 

.68985276+03 

.00000006+00 

.30 

.OOOOOOOE+OO 

.30000006-00 

. 7760843E+03 

. 1310063E- 10 

.40 

. OCOOOOOE+OO 

,40000006+00 

. 7884030E+C3 

.1007527E-10 

.50 

.OOOOCOOE+OO 

.50000006+00 

.78392356+03 

.OOOOOOOE+OO 

.60 

.00000006+00 

,60000006+00 

.77608436+03 

.OOOOOOOE+OO 

.70 

.OOOOOOCE+OG 

. 7000000E+00 

,763244?t+03 

.00000006+00 

,80 

.00000006+00 

.BDOOCOOE+OO 

.76121676+03 

.OGOOOOOE+OO 

.90 

.OOOOOCOE+OC 

. 9000000E+00 

.75510896+03 

.ocooooot+oo 

1.00 

.OOOOOOOE+OO 

.OOOOOOOE+OO 

.OOOOCOOE+OO 

.OOGOOOOE+OO 

AVERAGE 

OF THE FOCAL 

LENGTHS (meters) = 

= .677655E+01 

INCREMENTAL DISTORTION {%) RANGES FROM 
.. END THERMAL LENSING OUTPUT.. 

0.0 TO . 1 3 1 E 

-10 


86 


ORIGINAL PAGE IS 
OF POOR QUALITY 


BLOCK OF TEMPERATURE PROFILES FRCH INPUT FILE #2 


TEMPERATURE DISTRIBUTION 


RADIAL DIFFUSIVITY 2694706+02 

RADIUS OF R00 25QOOOE+GO 

LENGTH OF ROC 7600006+01 

TIME IN SECONDS ,3450136+01 


l/R 

.00 

.10 

.20 

.30 

.40 

.00 

.72332E+QQ 

. 72383E+0C 

. 72383E+00 

.723836+00 

.72383E+00 

.10 

.713306+00 

.71331E+Q0 

.71331E+00 

.713316+00 

.713316+00 

.20 

.694036+00 

.694 046+00 

.69404E+00 

.69404E+00 

.694046+00 

.30 

.663206-00 

.663206+00 

.663206-00 

.66320E-00 

.663206+00 

.40 

.62046E+00 

.620466+00 

.620466+00 

.62046E+C0 

.620466+00 

.50 

. 56521 E+00 

. 565226+00 

.565226+00 

.565226+00 

.565226+00 

.60 

.49683F+00 

.496S9E+0C 

.496696+00 

.496696+00 

.496896+00 

.70 

.415526+00 

.415526-00 

.41552E+00 

.415526+00 

.415526+00 

.60 

.32288E+00 

.322386+00 

.32286E+00 

.322866+00 

.322686+00 

.90 

.22190E+00 

.221906-00 

.221906+00 

.221906*00 

.22190E+00 

S.00 

. 1Q5566+CO 

.105566+00 

.105566+00 

.105566+00 

.105566+00 


THERMAL LENS! KG PROGRAM OUTPUT DATA 
INDEX OF REFRACTION MATRIX 


.00 

. 182506+01 

. 182506+01 

.182506+01 

.182506+01 

.132506+01 

.10 

182506+01 

. 1 82506 +C1 

.182506+01 

.182506+01 

. 18250E+01 

.20 

.182506+01 

.182506+01 

.132506+01 

.182506+01 

.182506+01 

.30 

.182506+01 

.182506-01 

.182506+01 

.182506+0 i 

.162506+01 

.40 

„ 132506+01 

.182506+01 

.182506+01 

.182506+01 

. 18250E+01 

,50 

.152506+01 

.182506+01 

.182506+01 

. 18250E+01 

. 18250E+01 

.60 

.182506+01 

.182506+01 

.182506+01 

.182506+01 . 

,132506+01 

.70 

.182506+01 

.182506-01 

. 1825CE+01 

. 18250E+01 

. 18250E+01 

,80 

. 18250E+01 

. 132506+01 

.182506+01 

.182506+01 

.182506+01 

.90 

.182506+01 

. 182506+01 

. 18250E+ni 

.182506+01 

.182506+01 

.00 

.182506+01 

.182506+01 

.182506+01 

.18250E+01 

.182506+01 


.50 

.60 

.70 

.80 

.90 

1.00 

.723836+00 

.723836+00 

.723836+00 

.723836+00 

.723836+00 

.723826+00 

.713316+00 

.713316+00 

.71331E+00 

.71 3316+00 

,713316+00 

.713306+00 

.694046+00 

.694046+00 

.694046+ 00 

.69404E+00 

.694046+00 

.694036+00 

.66320E+00 

.663206+00 

.663206+00 

.66320E+00 

.663206+00 

.663206+00 

.620476+00 

.62046E+00 

.620466+00 

.62046E+00 

.620466-00 

.620466+00 

.56522E+00 

.565226+00 

.565226+00 

.565226+00 

.565226+00 

.565216+00 

.496896+00 

. 496896+00 

.496896+00 

.49689E+00 

.496896+00 

.496886+00 

.415526+00 

.415526+00 

.415526+00 

.41552E+00 

.415526+00 

.415526+00 

.32288E+00 

.322886+00 

.322886+00 

.322886+00 

.322S8E+00 

.322886+00 

.221906+00 

.22190E+00 

.221906+00 

.221906+00 

.22190E+00 

.22190E+00 

.105566+00 

.105566+00 

.105566+00 

.105566+00 

.105566+00 

.105566+00 


182506+01 

,182506+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

182506+01 

.182506+01 

.182506+01 

. 18250E+D1 

. 18250E+01 

.182506+01 

182506+01 

.182506+01 

.162506+01 

. 18250E+01 

. 18250E+01 

.182506+01 

18250E+01 

.182506-01 

.182506+01 

. 1P.250E+01 

.13250E+01 

.182506+01 

18250E+01 

.182506+01 

.182506+01 

.182506+01 

.1S250E+01 

.182506+01 

18250E+01 

.182506+01 

. 18250E+01 

. 18250E+01 

. 18250E+01 

.182506+01 

1825QE+Q1 

.182506+01 

. 18250E+Q1 

.182506+01 

. 18250E+01 

.182506+01 

18250E+01 

.182506+01 

. 182506*01 

. 18250E+01 

.182506-01 

.182506+01 

182506+01 

.182506+01 

. 18250E+01 

.182506+01 

.182506+01 

. 18250E+01 

,182506+01 

I8250E+ 0 ! 

.182506+01 

.182506+01 

.182506+01 

. 18250E+01 

182506+01 

. 18250E+01 

.132506+01 

.182506+01 

.182506+01 

.162506+01 




THERMAL 

LENSING CAL CUAT IONS 




WAVEL6NG T H OF 

INTEREST (uu) 

. 1150C0E+01 



TEMPERATURE OF 
TIME (sec) 

i NT EREST (deg C> 

.375000E+02 

.345013E+01 


RAY 

ANGLE CUT 

RADIUS 

FOCAL 

DISTORTION 

- 

(no. ) 

(radians) 

(cn) 

LENG(cm) 

(X) 

/ 

.00 

.00000006+00 

.00000006+00 

.OOOOOOOE+OO 

. COOOOOOE+OO 


.10 

.OOOOOOOE+OO 

. 1000000E+00 

.34492636+03 

.57731566-11 


.20 

.00000006-00 

.20000006+00 

.6898527E+03 

. 1 393881 E- 10 


.30 

.00000006+00 

.zooQcooe+oc 

.7760843E+03 

.OQCOOOOE-OO 


.40 

.00000006+00 

.40000006+00 

.6898527E+03 

. 1333073E- 08 


.50 

.00000006+00 

.50000006+00 

.71859656+03 

* . 154765 iE- 10 


.60 

.00000006+00 

,50000006+00 

.73043236+03 

.22167456 -10 


.70 

.OOOOOOOE+OO 

.70000006+00 

.73484306+03 

.00000006+00 


.80 

.OCOOOOOE +00 

.80000006+00 

.7358428E+03 

.00000006+00 


.90 

.OOOOOOOE+OO 

.90000006+00 

.75510896+03 

OOOOOOOE+0.1 


1.C0 

.OOOOOOOE+OO 

.OOOOOOOE+OO 

.00000006*00 

.00C0G006+C0 


AVERAGE OF THE FOCAL LENGTHS (meters) = .68617SE+01 

INCREMENTAL DISTORTION (X) RANGES FROM 0.0 TO .1336-08 


.. END THERMAL L6NSING OUTPUT.. 


ORIGINAL PAGE IS 
OF POOR QUALITY 


BLOCK OF TEMPERATURE PROFILES FROM INPUT FILE #2 


TEMPERATURE DISTRIBUTION 


RADIAL D1FFUSIVITY 269470E+02 

RADIUS OF ROD 2500006+00 

LENGTH OF ROD 780000E+01 

TIME IN SECONDS 3500Q06+01 


L/R 

.00 

.10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

.90 

1.00 

.00 

.64627E+0Q 

.64627E+00 

.64627E+00 

.64627E+GQ 

.646276+00 

,646276+00 

.646276+00 

•64627E+00 

.646276+00 

.646276+00 

.646276+00 

.10 

.640066+00 

.64007E+0Q 

.640Q7E+00 

.640076+00 

.64007E+00 

. 640076+00 

.640076+00 

.64007E+00 

.640076+00 

.640076+00 

.640066+00 

.20 

.62144E+00 

.621456+00 

.621456+00 

.62145E+00 

.62145E+00 

.62145E+00 

.621456+00 

.62145E+00 

.62145E+00 

.62145E+00 

.621446+00 

.30 

.59041E+G0 

.59Q41E+00 

.590416+00 

.590416+00 

.59041E+00 

.59041E+00 

,590416+00 

. 5904 1 E+00 

.5904 IE +00 

.590416+00 

.59041 E+00 

.40 

.54693E+00 

.546946+00 

.546946+00 

.54694E+00 

.54694E+0G 

. 54694E+G0 

.54694E+00 

,546946+00 

.54694E+00 

.54694E+00 

.546936+00 

.50 

.49096E+00 

.490966+00 

.490966+00 

.490966+00 

.49096E+00 

. 49096E+00 

.49096E+00 

.490966+00 

.49096E+00 

.490966+00 

.490966+00 

.60 

.422446+00 

.42244E+00 

.42244E+00 

.422446+00 

,422446+00 

.42244E+00 

.422446+00 

.42244E+00 

.42244F+00 

.422446+00 

.422446+00 

.70 

.34169E+00 

. 341696+00 

. 34169E+00 

.341696+00 

.341696+00 

.34169E+00 

.34169E+00 

.341696+00 

.341696+00 

. 341696+00 

.34169E+00 

.80 

.25008E+00 

. 250086+00 

. 250086+00 

.25008E+00 

. 25008E+0Q 

.250086+00 

.250C8E+00 

.250086+00 

.250086+00 

.25008E+00 

.25008E+00 

.90 

. 15093E+00 

.150936+00 

. 15093E+00 

.150936+00 

♦15C93E+00 

. 1 5093E+OO 

. 150936+00 

.150936+00 

.15093E+00 

. 15093E+00 

. 15093 E+00 

1.00 

.49856E-01 

.49S56E-01 

.49856E-01 

.49857E-01 

.498576-01 

.49857E-01 

.49857E-01 

.498576-01 

.49856E-0I 

.4985 6E-01 

.49856E-01 


THERMAL LENS INS PROGRAM OUTPUT DATA 
INDEX OF REFRACTION HATRIX 


.00 

. 1025OE+O1 

. 13250E+01 

.182506+01 

.182506+01 

. 18250E+01 

. 18250E+01 

.182506+01 

.182506+01 

. 18250E+01 

. 18250E+01 

.182506+01 

.10 

. 18250E+01 

. 1 3250E+01 

. 18250E+01 

. 18250E+01 

. 18250E+01 

.182506+01 

. 18250E+01 

.182506+01 

.182506+01 

. 1B250E+01 

. 18250E+01 

.20 

. 18250E+01 

. 18250E+0 1 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

. 18250E+01 

.182506+01 

.132506+01 

. 18250E+01 

.18250E+01 

.30 

. 18250E+01 

.182506+01 

.182506+01 

. 182506+01 

. 1 8250E+01 

. 1 8250E+01 

. 18250E+01 

.182506+01 

.182506+01 

.182506+01 

. 18250E+01 

.40 

. 18250E+01 

.182506+01 

.182506+01 

. 18250E+01 

.182506+01 

. 1 8250E+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

.50 

. 18250E+01 

.182506+01 

. 18250E+01 

.182506+01 

.18250E+01 

. 18250E+01 

.182506+01 

. 18250E+01 

.182506+01 

.182506+01 

.182506+01 

.60 

. 182506+01 

. 18250E+C1 

. 1B2506+01 

. 18250E+01 

.182506+01 

. 18250E+01 

.182506+01 

,182506+01 

,182506+01 

.182506+01 

.182506+01 

.70 

, 18250E+01 

. 1 8250E+01 

.182506+01 

, 18250E+01 

.182506+01 

. 18250E + 01 

.182506+01 

. 18250E+01 

.182506+01 

. 18250E+01 

. 18250E+Q1 

.80 

. 18250E+01 

.182506+01 

.ie250£+01 

. 1825GE+01 

.182506+01 

.182506+01 

. 18250E+01 

. 1B250E+01 

.182506+01 

.182506+01 

. 16250E+Q1 

.90 

. 182506+01 

. 1 825C5+Q1 

.182506+01 

.182506+01 

. 18250E+01 

.182506+01 

. 18250E+01 

. 18250E+01 

. 18250E+01 

.182506+01 

. 18250E+01 

1.00 

. 18250E+01 

.182506+01 

.182506+01 

.162506+01 

. 18250E+01 

.182506+01 

. 18250E+01 

. 18250E+01 

.182506+01 

. 18250E+01 

. 18250E+01 


THERMAL LENSING CAtCUATIONS 


WAVELENGTH OF INTEREST(um) 1150006+01 

TEMPERATURE OF INTEREST(dcg C).. .3750006+02 

TIME (sec) 3500006+01 


RAY 

ANGLE CUT 

RADIUS 

FOCAL 

DISTORTION 

(no.) 

(rad? ans) 

(cm) 

LENG(cm) 

(X) 

.00 

.00000006+00 

.00000006+00 

.OOOOOOOE+OO 

.OOCOUOCC+OO 

.10 

.00000006+00 

.10000006+00 

.OOOOOOOE+OO 

.45519136-11 

.20 

.OOGOOOOE+OO 

. 2C000006+00 

.68985276+03 

-.12295726-10 

.30 

.OOOOOOOE+OO 

.3000000E+00 

.7760843E+03 

.00000006+00 

.40 

.00000006+00 

.40000006+00 

.7884030E+03 

.97977156-11 

.50 

.00000006+00 

.50000006+00 

.78392356+03 

.00000006+00 

.60 

.00000006+00 

.60000006+00 

.7760843E+03 

.00000006+00 

.70 

.00000006+00 

.70000006+00 

. 7682449E+03 

.00000006+00 

.80 

.00000006+00 

.80000006+00 

.7612167E+03 

.00000006+00 

.90 

.00000006+00 

.90000006+00 

. 7551089E+03 

.00000006+00 

.00 

.00000006+00 

.00000006+00 

.00000006+00 

.00000006+00 


AVERAGE OF THE FOCAL LENGTHS (meters) = .6776586+01 

INCREMENTAL DISTORTION (X) RANGES FROM 0.0 TO .9806-11 

.. END THERMAL LENSING OUTPUT.. 


ORIGINAL PAGE IS 
OF POOR QUALITY 


88 



BLOCK OF TEMPERATURE PROFILES FROM INPUT FILi M2 


TEMPERATURE DISTRIBUTION 


RADIAL DIFFUSIVTTY 2694706+02 

RADIUS OF ROO 2500006+00 

LENGTH OF ROO ......... 7600006+01 

TIME IN SECONDS 3500136+01 


L/r .00 .10 .20 .30 .40 .50 .60 .70 

.00 .723826+00 .723836+00 .723836-00 .723836+00 .72383E+00 .723836+00 .723836+00 .72383E+00 

* 10 .71330E+00 .713316+00 .713316+00 .71331E+00 .713316+00 .71331E-0G .713316+00 .71331E+00 

.20 .694036+00 .694046+00 .694046+00 .69404E+00 .694046+00 .694046+00 .69404E+00 .694046+00 

*30 ! 663206+00 .663216+00 .66321E+00 .663Z1E+00 .66321E+00 .66321E+00 -66321E+00 .683216+00 

.40 .620466+00 .620476+00 .620476+00 .620476+00 .62047E+00 .62047E+00 .62047E+00 .62047E+00 

’50 ’.56521E+00 .56522E>00 .56522E+00 .565226+00 .565226+00 .565226+00 .56522E+00 .565226+00 

.60 ,496886+00 .496896+00 .496896+00 .496896+00 .496896+00 .<<96896+00 .49689E+00 .496896+00 

*70 ^415526+00 .41552E+00 .415522*00 .415526+00 .415526+00 .415526+00 .415526+00 .415526+00 

*50 .322886+00 .322882*00 .32288E*00 .3228SE*00 .322886*00 .32288E*00 .322886+00 .322886+00 

.90 .22190E+00 .22190E+00 .221906*00 .22190E+00 .22190E+00 .22190E+00 .221906+00 .22190E+00 

1.00 . 10556E+00 .105566+00 .10556E+00 .105566+00 .10556E+00 .105566+00 .10556E+00 .105566+00 

THERMAL LENS INC PROGRAM OUTPUT DATA 
INDEX OF REFRACTION MATRIX 

.00 . 18250E+01 .1825CE+01 .18250E+01 .18250E+01 .T8250E+01 .182506+01 .182506+01 .182506+01 

,10 .182506+01 .182506+01 .18250E+01 .182506+01 .182506+01 .182506+01 .1025OE+O1 .18250E+01 

.20 . 18250E+01 , 18250E+01 .18250E+01 . 182506+0 1 .182506+01 .182506+01 .18250E+01 .182506+01 

.30 .1B250E+01 .182506+01 .18250E+01 .182506+01 .18250E+01 .182506+01 .18250E+01 .182506+01 

.40 .182506+01 . 18750E+01 . 18250E+01 .182506+01 .18250E+01 .182506+01 .18250E+01 .182506+01 

.50 ! 182506+01 . 18250E+01 .182506+01 .18250E+01 .18250E+01 .182506+01 .182506+01 .182506+01 

.60 .182506+01 .1825CE+Q1 ,182506+01 .18250E+01 .182506+01 .182506+01 .182506+01 .182506+01 

.70 .182506+01 .182506+01 .182506+01 .182506+01 .182506+01 .1B250E+01 .182506+01 .182506+01 

.80 .182506+01 .182506+01 .1B250E+01 .18250E+01 .18250E+01 .18250E+01 .182506+01 .18250E+01 

.90 . 18250E+01 .182506+01 .182506+01 .182506+01 .182506+01 .18250E+01 .182506+01 .182506+01 

1 .00 .182506+01 . 182506 *-0 1 .182506+01 .182506+01 .182506+01 .182506+01 ,182506+01 .182506+01 


.80 

.90 

1.00 

.72383E+00 

.723836+00 

.72382E+00 

.713316+00 

.713316+00 

.71330E+00 

.694046+00 

.694046+00 

.69403E+00 

.663216+00 

.663216+00 

.663206+00 

.620476+00 

.62047E+00 

.620466+00 

.565226+00 

.56522E+0G 

.565216+00 

.496896+00 

.496896+00 

.496886+00 

.41552E+00 

.415526+00 

.41552E+00 

.32288E+00 

.32288E+00 

.322886+00 

.221906+00 

.221906+00 

.22190E+00 

.105566+00 

.105566+00 

.105566+00 


182506+01 

. 1 8250E + 01 

.182506+01 

182506+01 

. 18250E+01 

.182506+01 

182506+01 

.182506+01 

.182506+01 

182506+01 

.182506+01 

.182506+01 

18250E+01 

.182506+01 

. 18250E+01 

18250E+01 

.182506+01 

.182506+01 

18250E+01 

.182506+01 

.182506+01 

182506+01 

. 18250E+01 

.182506+01 

13250E+01 

.182506+01 

.182506+01 

18250E+01 

.182506+01 

.162506+01 

182506+01 

.182506+01 

.182506+01 


THERMAL LENS! NG CALCUATIONS 



WAVELENGTH OF 

INTEREST (um) . . 

1150006 

+C1 


TEMPERATURE OF 

I NTEREST (deg 

C>.. .375000E+02 


TIME (sec) 


350013E+01 

RAY 

ANGLE OUT 

RADIUS 

FOCAL 

DISTORTION 

(no.) 

(radians) 

(cm) 

LENG(cm) 

(X; 

.00 

.00000006+00 

.00000006+00 

.OOOOOOOE+OO 

.00000006+00 

.10 

.00000006+00 

.10000006+00 

.34492636+03 

.57731566-11 

.20 

.00000006+00 

. 2000000E+QU 

.68985276+03 

- .47134486- 1 1 

.30 

.00000006+00 

.30000006+00 

.77608436+03 

. 355271 2E-11 

.40 

.00000006+00 

. 4 000000 E+ 00 

.6898527E+03 

.9936493E-11 

.50 

.00000006+00 

.5000000E+00 

.71859656+03 

-.66613386-12 

.60 

.OOOOOOQE+OQ 

.60000006+00 

.73C4323E+03 

.2216745E-10 

.70 

.00000006+00 

.70000006+00 

.73484306+03 

.00000006+00 

.80 

.OOCOOOOE+DO 

.80000006+00 

, 7358428E*03 

.OOOOOOOE+OO 

.90 

.QOOQOOOE+OO 

.90000006+00 

.75510896+03 

.00000006+00 

.00 

.OOOOOOOE+OO 

.OOOOOOOE+OO 

.OOOOOOOE+OO 

.00000006+00 


AVERAGE OF THE FOCAL LENGTHS (meters) = .6861716+01 

INCREMENTAL DISTORTION (X) RANGES FROM 0.0 TO .2226*10 


.. END THERMAu LEN31NG OUTPUT . . 
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BLOCK OF TEMPERATURE PROFILES FROM INPUT FILE #2 


TEHPERATURE DISTRIBUTION 

RADIAL D I FFUSIVITY 269470E+02 

RADIUS OF ROO 2500006+00 

LENGTH OF ROO 760000E+01 

TIME IN SECONDS 3550006+01 


l/N 

.00 

.10 

.20 

.30 

.40 

.00 

.64627E+00 

. 64627E+00 

.64627E+00 

.64627E+00 

.646286+00 

.10 

.640066+00 

.640076+00 

. 64 0G7E+0Q 

.64Q07E+00 

.64Q07E+QQ 

.20 

.621446+00 

.62145E+00 

.62145E+00 

.621456+00 

.62145E+D0 

.30 

.59041E+00 

.59041E+00 

.59041E+00 

.59041E+00 

.59042E+00 

.40 

.546946+00 

.546946+00 

.54694E+00 

.54694E+00 

. 54694E+00 

.50 

.490966+00 

.49097E+00 

.49097E+00 

.49097E+00 

.49097E+00 

.60 

.42244E+00 

.42244E+0G 

•42244E+00 

.4Z244E+00 

. 42244E+00 

.70 

.34 1696+ DO 

.34169E+00 

.341696+00 

.341696+00 

.34 169E+0Q 

.80 

.250086+00 

.25Q08E+Q0 

.250086+00 

. 250086+00 

. 25008E+00 

.90 

. 15093E+00 

. 15093E+00 

. 15093E+00 

.150936+00 

. 15093E+Q0 

1.00 

.49856E- 01 

.498575*01 

.49B57E-01 

.498575-01 

.49857E-01 


THERMAL LENSING 

FROGRAM OUTPUT DATA 




INDEX OF REFRACTION MATRIX 


.00 

. 1825OE+01 

. 1825GE+01 

.182506+01 

. 18250E+01 

.182506+01 

.10 

. 18250E+01 

. 18250E+01 

. 182506+01 

. 1825QE+01 

. 182506+01 

.20 

. 18250E+01 

. 18250E+01 

. 18250E+01 

. 18250E+01 

. 18250E+Q1 

.30 

. 18250E+01 

. 1B250E+Q1 

. 182506+01 

. 18250E+Q1 

. 18250E+01 

.40 

. 1825GE+01 

.18250E+01 

. 18250E+01 

. 182506+01 

. 18250E+01 

.50 

. 182506+01 

. 18250E+G1 

.182506+01 

. 18250E+01 

. 182506+01 

.60 

. 182506+01 

. 1S250E+C1 

. 182506+01 

. 18250E+01 

.182506+01 

.70 

. 182506+01 

. 18250E+01 

. 18250E+01 

. 18250E+O1 

. 18250E+01 

.80 

.182506+01 

. 1 82506+01 

. 18250E+01 

. 182506+01 

. 18250E+01 

.90 

. 1825DE+01 

. 182506+01 

.182506+01 

. 18250E+Q1 

, 1825QE+01 

1.00 

.162506+01 

. 1825Q6+0T 

. 18250E+G1 

. 18250E+Q1 

. 1025OE+Q1 


.50 

.60 

.70 

.80 

.90 

1.00 

.64628E+00 

.6462BE+00 

.646275+00 

•64627E+00 

.64627E+00 

.64627E+00 

.64007E+00 

.64007E+00 

. 64 00 7E+ 00 

.64007E+00 

.64007E+00 

.64006E+00 

.62145E+00 

.62145E+Q0 

.62145E+00 

.62145E+00 

.621456+00 

.62144E+00 

.590426+00 

.59042E+00 

.59041E+00 

■59041E+00 

.59041E+00 

-59041E+00 

.54694E+00 

.54694E+00 

.54694E+00 

.54694E+00 

.54694E+00 

. 54694 E+00 

. 49097E+00 

.49097E+00 

.49097E+00 

.49097E+00 

. 49097E+00 

.49096E+00 

.42244E+00 

.42244E+00 

•42244E+00 

. 42244E+00 

.42244E+00 

.42244E+00 

.34169E+00 

.34169E+00 

.34169E+00 

.34 169E+00 

. 34 169E+00 

.34169E+00 

. 250086+00 

.25008E+00 

.25008E+00 

. 25008E+00 

. 2500SE+00 

.25008E+00 

. 15093E+00 

. 15093E+00 

. 15093E+00 

.150935+00 

.15093E+00 

.15093E+00 

49857E-01 

.49B57E-01 

.49857E-01 

-49857E-01 

.49857E-01 

.498566-01 


.182506+01 

.182506+01 

.182506+01 

. 18250E+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

.18250E+01 

.182506+01 

.182506+01 

.182506+01 

. 18250E+01 

. 182506+01 

. 182506+01 

.182506+01 

.182506+01 

•18250E+01 

,182506+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

.18250E+01 

,182506+01 

.182506+01 

.182506+01 

. 1S250E+01 

.182506+01 

. 18250E+01 

,182506+01 

.182506+01 

.182506+01 

. 1 82506+01 

.182506+01 

.182506+01 

182506+01 

. 18250E+01 

.182506+01 

.182506+01 

, 18250E+01 

. 18250E+01 

182506+01 

.182506+01 

.152506+01 

. 1825QE+C1 

.182506+01 

.182506+01 

18250E+01 

. 18250E+01 

. 182506+01 

.182506+01 

. 18250E+01 

.18250E+01 

182506+01 

.182506+01 

.182506+01 

. 18250E+01 

.182506+01 

.182506+01 

182506+01 

.182506+01 

. 18250E+01 

.182506+01 

. 18250E+01 

. 182506+01 


THERMAL LENS I HG CALCUATIONS 

WAVELENGTH OF I NTERESTCun) . .1150006+01 

TEMPERATURE OF INTEREST (deg C).. .375000E+Q2 

TIHE (sec) 3550006+01 


RAY 

ANGLE OUT 

RADIUS 

FOCAL 

DISTORTION 

{no. ) 

(radians) 

(cm) 

LENG(cm) 

(X) 

.00 

.00000006+00 

.OOOOOOOE+OO 

.OOOOOOOE+OO 

.00000006+00 

.10 

.00000006+00 

. 1 OOOCOOE+OO 

.OOOOOOOE+OO 

.45519136-11 

.20 

. QQQOOOOE+QU 

. 2000000E+00 

.6893527E+03 

*.12295726-10 

.30 

. OOQOOOOE+OO 

.3000000E+00 

. 7760843E+03 

.99635116-09 

.40 

. OOOOOOOE+OO 

. 4000000E+00 

, 7884 0306+03 

.OCOOOQQE+OO 

.50 

. DOOOOOOE+OO 

. 500Q000E+C0 

.78392356+03 

- . 1709743E-1 1 

.60 

.OOOOOOOE+OO 

.6000000E+00 

. 7760843E+03 

.00000006+00 

.70 

.00000006+00 

.70000006+00 

.76824496+03 

.00000006+00 

.80 

.00000006+00 

.8000000E+00 

.7612167E+03 

. COQQOOOE+OO 

.90 

.00000006+00 

.9000000E+00 

.75510896+03 

.00000006+00 

1.00 

.OOOOOOOE+OO 

.00000006+00 

.00000006+00 

.OOOOOOOE+OO 


AVERAGE OF THE FOCAL LENGTHS (meters) = .677658E+01 

INCREMENTAL DISTORTION (X) RANGES FROM 0.0 TO .9966-09 

. . END THERMAL LENSING OUTPUT. . 
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BLOCK OF TEMPERATURE PROFILES FROM INPUT FILE #2 


TEMPERATURE DISTRIBUTION 


RADIAL D I FFUSIVITY 269470E+02 

RADIUS OF ROC 250000E+G0 

LENGTH OF ROO 760000E+01 

TIME IN SECONDS 3550136+01 


L/R 

.00 

.10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

.90 

1 .00 

.00 

. 72382E+00 

. 723836+00 

.72383E+00 

.72383E+00 

.72383E+00 

.72383E+D0 

.72383E+00 

.72383E+00 

.723836+00 

.723836+00 

.723826+00 

.10 

.713316+00 

.713316+00 

.71331E+00 

.71331E+00 

.71331E+00 

.713316+00 

.71331E+00 

.713316+00 

.71331E+00 

.713316+00 

.71331E+00 

.20 

.694036+00 

. 69404 E+ 00 

. 69404E+00 

.694046+00 

.69404E+00 

.69404E+00 

.69404E+0Q 

.69404E+00 

.69404E+00 

.694046+00 

.69403E+00 

.30 

.66320E+00 

.66321E+00 

. 66321 E +0C 

.663216+00 

.66321E+00 

. 66321E+00 

.663216+00 

.66321E+00 

.66321E+00 

.663216+00 

.663206+00 

.40 

.620466+00 

. 62047E+C0 

.620476+00 

.620476+00 

.62047E+00 

.62047E+00 

.62047E+00 

.620476+00 

.62047E+00 

.6204 7E+00 

.62046E+00 

.50 

.565216+00 

.365226+00 

. 56522E+00 

.565226+00 

.56522E+00 

. 56522E+00 

.565226+00 

.565226+00 

. 56522E+00 

.565226+00 

.56521E+00 

.60 

.49688E+00 

.496896+00 

.49689E+00 

.496896+00 

.49689E+00 

.496896+00 

.496896+00 

.496896+00 

.496896+00 

.49689E+00 

.496896+00 

.70 

.415526+00 

.41552E+00 

.41552E+00 

.415526+00 

.4 1552E+-00 

.415526+00 

.415526+00 

.415526+00 

.4 1552E+00 

.415526+00 

.415526+00 

.80 

.322886+00 

.32288E+00 

.322886+00 

.322886+00 

.322886+00 

.322886+00 

.322886+00 

.322586+00 

.322886+00 

.322586+00 

.322886+00 

.90 

.22190E+00 

.22190E+00 

.221906+00 

.22190E+00 

.221906+00 

.221906+00 

.221906+00 

.221906+00 

.221906+00 

.221906+00 

.22190E+00 

.00 

. 1C556E+00 

.10556E+00 

. 10556E+0C 

.10556E +00 

.105566+00 

.105566+00 

.105566+00 

. 105566+00 

.105566+00 

.105566+00 

.105566+00 


THERMAL LENS INC PROGRAM OUTPUT DATA 


INDEX OF REFRACTION MATRIX 


.00 

.182506+01 

. 18250E+01 

.132506+01 

. 18250E+01 

. 18250E+01 

. 18250E+01 

.182506+01 

. 18250E+01 

.182506+01 

.182506+01 

.182506+01 

.10 

.182506+01 

.182506+01 

.182506+01 

.18250E+01 

.182506+01 

.182506+01 

.182506+01 

. 18250E+01 

.182506+01 

.182506+01 

. 18250E+C1 

.20 

.132506+01 

.182506+01 

.18250E+01 

.182506+01 

.182506+01 

. 18250E+01 

.182506+01 

. 18250E+01 

.182506+01 

.182506+01 

.182506+01 

.30 

.182506+01 

. 182506 ■‘■01 

.182506+01 

.182506+01 

. 18250E+01 

. 18250E+01 

.182506+01 

.18250E+01 

.182506+01 

. 1825C6+01 

.182506+01 

.40 

.182506+01 

. 1825OE+01 

. 18250E+01 

.18250E+01 

.18250E+01 

.182506+01 

.182506+01 

.182506+01 

. 18250E+01 

.182506+01 

.18.2506+01 

.50 

. 18250E+01 

.182506+01 

.182506+01 

. 18250E+01 

. 1 S250E+01 

.182506+01 

.182506+01 

.182506+01 

. 1S250E+01 

. 18250E+01 

. 18250E+01 

.60 

.182506+01 

.182506+01 

. 13250E+01 

.182506+01 

.182506+01 

.182506+01 

. 18250E+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

,70 

.182506+01 

.182506+01 

. 1B250E+ 01 

. 18250E+01 

. 18250E+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

. 1825G6+01 

. 13250E+01 

.80 

. 1S250E+ 01 

.182506+01 

. 18250E+01 

.18250E+01 

. 18250E+01 

.182506+01 

.182506+01 

. 18250E+01 

. 18250E+01 

. 182506+01 

. 18250F+01 

,90 

.182506+01 

.182506+01 

. 18250E+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

. 182506*01 

.182506+01 

.182506+01 

,00 

.182506+01 

. 182506+01 

.1B250E+01 

.182506+01 

.182506+01 

.182506+01 

.182506+01 

. 1B250E+01 

.182506*01 

. 18250E+D1 

.182506+01 


THERMAL LEWS 1 WG CALCUATIONS 



WAVELENGTH OF 

INTEREST(um) 

.1150006+01 



TEMPERATURE OF 

I NTERFST (deg C). 

.3750006+02 



TIME (see) 


.355013E+01 


RAT 

ANGLE OUT 

RADIUS 

fOCAL 

DISTORTION 

{no. 

) (radians) 

(cm) 

LENG(cm) 

(50 

.00 

.00000006+00 

.00000006+00 

.00000006+00 

.00000006+00 

.10 

.00000006+00 

.10000006+00 

.34492636+03 

.00000006+00 

.20 

.oooooooc+cc 

.2000000E+00 

.68985276+03 

. 25535 1 6E - 1 1 

.30 

.00000006+00 

.30000006+00 

.77608436+03 

,35527126-11 

.40 

.OOOOOOOE+GO 

.40000006+00 

.68935276+03 

.99364936-11 

.50 

.OOCOOOOE+CO 

.50000006+00 

.71859656+03 - 

.66613386-12 

.60 

.00000006+00 

.60000006+00 

.7304323E+C3 

2216745E* 10 

.70 

.00000006+00 

.70000006*00 

.7348430E+G3 

OOOOOOOE+OO 

.80 

•OOOOOOOE+OO 

.80000006+00 

.735B428E+C3 

OOOOOOOE+OO 

.90 

.OOOOOOOE+CO 

.90000006+00 

. 7551039E+03 .OOOOOOOE+OO 

1.00 

.00000006+00 

.00000006+00 

.OOOOOOOE+OO 

OOOOOOOE+OO 


AVERAGE OF THE FOCAL LENGTHS (meters) = .656171E+01 

INCREMENTAL DISTORTION (X) RANGES FROM 0.0 TO .222E-10 

.. END THERMAL LEWSING OUTPUT.. 


ORIGINAL PAGE IS 
CF POOR QUALITY 


91 



References 


1. Farrukh, U. 0., Buoncristiani , A. M. , and Byvik, C. E., 
An Analysis of the Temperature Distribution in Finite 
snl id-State Laser Rods , IEEE Journal of Quantum 
Electronics, Vol . 24, No. 11, November 1988. 

2. Koechner , Walter., Solid State Laser — Eng i n e e r i n£L , 

Springer- Verlag : Berlin, Heidalburg, New York, 1988. 

3. Farrukh, U. 0., Buoncristiani, A. M. , and Byvik, C. E., 
Time Dependent Temperature Distri bution in Pulsed Ti : 
Sapphire Lasers , Proceedings of SPIE-The International 
Society for Optical Engineering, Vol. 892, 1988. 

4. Bendow, B. and Gianino, P. D., Thermal — Lensing of — Laser 
Beams in Optically Transmitting Mater ials - I. .. General 
Formul ation , Journal of Applied Physics, Vol. 2, pp. 
1-10, 1973. 

5. Koechner, W., Transient Thermal Profiles — in Op _t ically . 
Pumped Laser Rods , Journal of Applied Physics, Vol. 44, 
No. 7, pp. 3162-3170, 1973. 

6. Antipenko, B. M. , Cross-Relaxation Schemes — for — Pumpin g. 
Laser Rod s, Soviet Physics, Technical Physics, Vol. 29, 
No. 2, pp . 228-230, 1984. 

7. Farrukh, U., Brackett, V., and Brockman, P., Steady- 
State Temperature in End-Pumped Repetitive ly Pulsed 
Laser Rods, Technical Digest, Paper presented at the Optical 
Society of America, 1989 Annual Meeting, 15-20 October 1989. 


8. Foster, J. D. and Osterink, L. M. , Thermal Effects in 
Nd:YAG Lasers , Journal of Applied Physics, Vol. 41, No. 
9, pp. 3656-3663, 1970. 

9. Koechner, W., Thermal Lensing in N d:Y AG _Las er Ro d,/ 
Journal of Applied Optics, Vol. 9, No. 11, PP- 2548 
2553, Nov. 1970. 

10. Sparks, M., Optical Distortion by H eated Windows — in. 
H igh Power Laser Systems . Journal of Applied Physics, 
Vol. 42, No. 12, pp. 5029-5046, 1971. 

11. Jasperse, J. R. and Gianino, P. D. , Th ermal Lensing in 
Infrared T.aser Window Materials , Journal of Applied 
Physics, Vol. 43, No. 4, pp . 1686-1693, 1972. 


92 



References (cont.) 


12. Bendow, B. and Gianino, P. D., Thermal Lensinq of Laser 

Beams in Optica l ly Transmitting Materials - II. Numer- 
ical Computation s, Journal of Applied Physics, Vol . 2, 

pp. 71-90, 1973. 

13. Carslaro, H. S. and Jaeger, J. C., Conduction of Heat 
in Solids , Oxford University Press: London, 1948. 

14. Kogelnik, H., Imaging of Optical Modes - Resonators 
with Internal Lenses , Bell Systems Technical Journal, 
Vol. 44, pp. 455-494, 1965. 

15. Timoshenko, S. and Goodier, J. N., Theory of 
Elasticity , McGraw-Hill: New York, 1951. 

16. Born, M. and Wolf, E., Principles of Opti cs, Pergamon 
Press: London, 1965. 

17. Koechner, W. and Rice, D. K., Effect of Birefrin g ence 
on the Performance of Linearly Polarized YAG: Nd 
Lasers , IEEE Journal of Quantum Electronics, Vol. QE6 , 
No. 9, pp . 557-566, September 1970. 

18. Koechner, W., Absorbed Pump Power, Thermal. Profiles and 
Stresses i n a CW Pump ed Nd : YAG Crystal , Journal of 
Applied Optics, Vol. 9, No. 6, pp . 1429-1434, June 
1970 . 

19. Yariv, A. and Yeh, P., Optical Waves in Crystals , John 
Wiley & Sons: New York, 1984. 

20. Wi 1 1 iams-Burg , J. A. and Barnes, N. P., Laser Perfor- 
mance, Thermal Focusing and D epolariza ti on Effec ts i n 

Nd:Cr:CSGG and Nd:YAG , Presented at SPIE Conference in 
Los Angeles, Ca . , January 1990. 


93 


NASA 

NWOn.ii At^Toaok s 
l' AOPrf*sirdf*jn 


1. Report No 

NASA CR-4300 


4. Title and Subtitle 


Report Documentation Page 


2. Government Accession No. 


Modeling of Thermal Lensing in Side and End-Pumped 
Finite Solid-State Laser Rods 


3. Recipient's Catalog No. 


5. Report Date 

June 1990 

6. Performing Organization Code 


7. Author(s) 

Vincent G. Brackett 


9. Performing Organization Name and Address 

Hampton University 
Physics Department 
Hampton, VA 23668 


12. Sponsoring Agency Name and Address 

NASA Langley Research Center 
Hampton, VA 23665-5225 


8. Performing Organization Report No 


10. Work Unit No. 

590-31-31 


11. Contract or Grant No. 

NAGl-666 


13. Type of Report and Period Covered 

Contractor Report 

14. Sponsoring Agency Code 


IS. Supplementary Notes 

Langley Technical Monitor: Philip Brockman 

Thesis presented to The Graduate College, Hampton University, in partial 
fulfillment of the requirements for the Degree of Master of Science. 


16. Abstract 


An analytical expression for approximating the time-dependent thermal 
focal length In finite solid-state laser rods has been derived. The 
analysis Is based on the temperature variation of the material refractive 
index caused by optical pumping of these rods. Several quantities were 
found to be relevant to this analysis. These quantities were the specific 
thermal profiles of the rods, type of optical pumping employed, type of 
cooling scheme employed (side and end-cooling parameters), and the 
specific material characteristics of the rods. The Thermal Lensing Model 
was formulated using the geometric ray tracing approach. The focal 
lengths are then approximated, by calculating the phase shift in the Index 
of refraction, as the different rays of an Incident plane wave are tracked 
through a lens-like crystal medium. The approach also applies In the case 
of Gaussian or parabolic pump beams. It Is shown, that the prediction of 
thermal focal length is in good quantitative agreement with experimentally 
obtained data. 


17. Key Words (Suggested by Author! s)) 

optical pumping, solid-state lasers, 
optical properties, rods, temperature 
effects, pulsed lasers 


18. Distribution Statement 

Unclassified - Unlimited 
Subject Category 36 


19. Security Classif. (of this report) 

Unclassified 


20. Security Classif. (of this page) 

Unclassified 


I 21. No. of pages 


22. Price 


NASA FORM 1626 OCT W 

For sale by the National Technical Information Service, Springfield, Virginia 22161-2171 


NASA- Langley, 1990 



